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The dynamical properties of hole motion in an antiferromagnetic background are determined in 
. . . one dimensional models in zero magnetic field, where spin isotropy holds, as well as in an external 

magnetic field. The latter case is also relevant, via particle-hole transformation, to the problem of 
hole propagation in one dimensional "superconductors". The singularities in the spectral function 
■ are investigated by means of bosonization techniques and perturbation theories. Results are then 

compared with Bethe ansatz solutions and Lanczos diagonalizations. The formalism also leads to 
^Li, interesting connections to the single impurity problem in Luttinger liquids. A rich structure is 

found in the spectral function whenever spin isotropy is broken, suggesting the presence of exotic 
momentum dependence in photoemission spectra of (quasi) one dimensional materials. 

(N ' 
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^ ■ I. INTRODUCTION 

, The spectral properties of a single hole in a quantum antiferromagnet still represent an outstanding problem in 

' the physics of strongly correlated electron systems. Although more than 25 years elapsed since the seminal work 

f — ■ ' of Brinkman and Rice ^ (BR) there is still no consensus on the nature of hole motion (coherent or incoherent) or 

On on the features of the long range distortion induced by the hole on the antiferromagnetic ordering. On the other 
hand, a full understanding of the dynamics of a single hole is clearly required before the problem of the hole-hole 

^ I effective interaction, mediated by the magnetic background, can be addressed. This issue, relevant in the low doping 

y ' regime, is a key problem in the framework of high temperature superconductivity. Furthermore, recent developments 
in angle resolved photoemission and inverse photoemission experiments |^-^, have made possible to extract the 

^ ' momentum dependent spectral function in several compounds, including high temperature superconducting materials 

Q I at stoichiometric composition, which are good quantum antiferromagnets. The photoemitted electron leaves a mobile 

O ' hole in the spin background: therefore these experimental studies directly address the problem of hole propagation in 

^ systems where electron correlations play a key role. 

• I— I , A widely accepted model to describe the physics of a quantum antiferromagnet is the well known Hubbard model 

■ 1^ at half filling (one electron per site): 

S ■ 

H=-tJ2(4.C,a+h.C.) + UJ2i^,^~l)in^l^l) (1-1) 

<i.j> i 

where {cia) is the creation (annihilation) operator of an electron with spin a at the lattice site i and the symbol 
< i,i > indicates nearest neighbor summations over an hypercubic bipartite lattice in arbitrary spatial dimension d. 
Henceforth even (odd) values for i indicate conventionally one of the two sublattices. The operator riia- = c\^Cia is 
the number operator of a particle with spin a at the given site i. This model is defined in a finite lattice with L sites 
and standard periodic boundary conditions. When the total number of particles = -t- equals the number of 
sites, (i.e. at p = Nc/ L = 1) this Hamiltonian is believed to develop a gap in the charge excitation spectrum. This is 
actually rigorous in one dimension where the exact Lieb and Wu solution [pi yields a finite gap for arbitrary repulsion 

2L 

4 ^ (-l)^c.T (1.2) 



U > and magnetization per site /i = — Vr— ^- A particle hole transformation 



maps the half filled Hubbard model at [/ > and magnetization /i into the same model at ^' — 0, negative interaction 
([/' = —U) and density p' = 1 — 2/i. This mapping also shows that the problem of a single hole in an antiferromagnet 
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at non zero magnetic field is relevant to understand the photoemission spectra in superconductors: in fact, there is 
quite a robust numerical evidence [Q that the negative U model is an s-wave superconductor in d = 2, and in one 
dimension the exact Bethe ansatz solution predicts quasi-long range order in the ground state ^ . 

Although the treatment of a single hole {Nc — L — 1) might seem a major simplification, there are only few results 
valid and accepted in more than one dimension. The Nagaoka theorem is a remarkable exception, stating that for 
d > 1 the ferromagnetic state with maximum total spin S* = ^(L — 1) is the unique ground state of the infinite U 
Hubbard model, apart for the trivial degeneracy of the 25* + 1 spin components 

At strong coupling the Hubbard model at half filling is mapped into a standard Heisenberg model with antiferro- 
magnetic superexchange J = The presence of a single hole modifies slightly this mapping: each site is singly 
occupied and the hole hops from site to site weakening local antiferromagnetic correlations. This process is described 
by the so called t — J model: 

H ^-t {c^Cja + h.c.) + J 'Y {S^■ Sj ~ ^fiinj). (1.3) 

<ij>,cr <i,j> 

where the constraint of no double occupancy is understood. The exact mapping from Hubbard to t — J also includes 
a three site term, which is neglected here because is believed not to change the physics of the model, at least in 
one dimension pO| , pl| . The last, density dependent, contribution in Eq. (1.3) can be also dropped for our purposes 



because it is effective only when more than one hole is present in the system. 

This work extends and develops the analysis presented in a previous Letter about the spectral properties of hole 
motion in one dimensional models of correlated electrons. In such a case, spin charge decoupling allows to describe 
the low energy physics by an effective Hamiltonian H written as the sum of two commuting parts H = + Hp, the 
former governing the spin degrees of freedom (spinons) and the latter the charge ones (holons). 

This picture is by no means new in the field of one dimensional electron systems: standard analytical treatments 
show that, if the excitation spectrum of the model is gapless, the low energy physics in both charge and spin sectors 
is described by a Luttinger liquid model [l"^-|l5|] for generic microscopic Hamiltonians. The case is different in the 
Hubbard model for positive U at half filling or, for negative U, at zero magnetic field and arbitrary density: Only 
one of the two sectors is gapless and the renormalization group (RG) equations of the generic (the so called g— ology) 
model flow to strong coupling. The g— ology model parameterizes the most general low energy interaction present in 
one dimensional, translationally invariant systems of spin one-half electrons. The model depends on several coupling 
constants gi, with i — 1, . . .4, which may also have a spin dependence || and _L for electrons interacting with the 
same or with opposite spins, respectively JTst . Conventionally, gi refers to backward scattering, 172 f^nd g4 to forward 
scattering and 173 to Umklapp scattering, the latter present only in lattice models at commensurate fillings. When the 
coupling constants g's go to strong coupling under the RG flow, they should cross the exactly solvable Luther-Emery 
line leading to a spin gap for g^ < \gi±\ and a charge gap for —2g2 < {gsl- The Hubbard model at half filling 
belongs to the first class at negative U and to the second at positive U, thereby providing a model Hamiltonian which 
encompasses the most general strong coupling fixed points. Therefore, understanding the Hubbard model at half 
filling and arbitrary magnetization would shed light on the physical behavior at the Luther-Emery fixed point. This 
would be particularly valuable for the dynamical properties because the retarded Green function 

G^ip,Lj)^- <-9\clAi^-H + Eo + ivr'cp^\'f> (1.4) 

of the Luther-Emery model is not exactly known, although a widespread prejudice ascribes no interesting features to 
this correlation function |^,|l^ due to the presence of a gap in the excitation spectrum. 

A first study of the single hole Green function in quantum antiferromagnets (QAF) was performed by Brinkman 
and Rice Q who calculated the Green function of the t — J model by neglecting quantum fiuctuations, i.e. by replacing 
the spin exchange interaction by its Ising form: Si ■ Sj SfS^ in Eq.(p_.3|). In the absence of holes, the Neel state is 
the classical ground state of the model, whereas the lowest spin excitation has a gap ~ J. This Hamiltonian can be 
thought of as the strong coupling limit of a model characterized by a gap both in the spin and in the charge sectors 
(say the half filled Hubbard model with uniaxial spin anisotropy). The problem of hole motion in this system has 
been studied in any dimensions, but its solution in d = 1 is particularly simple and instructive. When the hole hops 
in the lattice it leaves a defect in the Neel background. As a result, the spectral weight A{ijj,p) — i ImG(p, w) shows 
a delta function contribution at the lowest excitation energy with non zero quasiparticle weight Z, together with an 
incoherent band separated by a gap. A brief discussion of these results is contained in Section |ll[ The situation is 
similar in the limit of infinite dimensions where there is no incoherent contribution and only a series of 5 function 
peaks at higher energies is left. The inclusion of quantum fluctuations, however, drastically changes this simple picture 
due to the presence of gapless excitations in the magnetic background. This is the subject of the present study which 
we organized as follows. 
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In Section II] an extremely useful mapping between the one hole Hamiltonian and an effective spin problem is given 
in some detail. In Section [V we discuss an interesting relationship, which emerges from the previous formulation, 
between hole dynamics and the impurity problem in Luttinger liquids. The main result of the present work can be 
summarized in the general structure we find in the Green function of all the models we have examined: 



G{p,t) = J ^Gh{p + Q,T)Zp{Q, 



-) 



(1.5) 



where Gh(k,t) is just a free propagator for the holon: lmGh{k,uj) = TrS{uj — eh{k)), eh being the holon dispersion 
energy. The function Z{Q,t) is completely determined by the spinon gapless excitations, and is highly non trivial 
with momentum dependent power law singularities and branch cuts. This decomposition is introduced in Section 
^ which also contains a discussion of the analytical properties of the function Z(Q, t). F ew specific examples, i.e. 
t — JxY model and the Bethe ansatz soluble models, are presented in Sections ^ and [VIl| respectively. A num erical 
evaluation of the non universal features of the spectral function in the Hubbard model is carried out in Section VIII 
while some conclusions are drawn in Section IX. 



II. HOLE MOTION IN ISING ANTIFERROMAGNETS 



The problem of hole motion in an Ising antiferromagnet can be solved exactly and provides the simplest model of 
hole dynamics in a magnetic environment characterized by a discrete excitation spectrum. This feature turns out to 



be responsible for the different behavior between the Ising and the XY case which will be discussed in Section [V] . 
The Hamiltonian describing hole hopping in an Ising model is a simple generalization of the t — J model where only 
the component of the spins is retained Q : 



H 



E 



[4' 



(2.1) 



with J > 0. We consider a chain with an even number L of sites and periodic boundary conditions. At half filling 
the ground state is a classical Neel state \J\f > . Out of the two states obtained by translation of one lattice vector 
we select the one with a spin down at the origin. The annihilation of the spin down electron at the origin defines our 
starting state |0 >= cqJA/" >. The Hamiltonian ( p.l[ ) acting on |0 > generates states which, in the thermodynamic 
limit, can be uniquely labeled by the position of the hole \R >. In fact, for L ^ oo, closed paths of the hole along 
the ring can be neglected and the retraceable path approximation becomes exact pj . Then, it is easy to check that 



H\R>=-t\\R-l> +\R+1 >] 



\R> ~'^5R,f,\R> 



(2.2) 



where Ej^ is the energy of the Neel state. The corresponding eigenvalue equation can be easily solved in this subspace: 
In the thermodynamic Hmit, the energy spectrum consists of a bound state with energy Ei, = Ej^ + ^^J — + l&t"^ 

and exponentially localized wavefunction 4'bi'r)- Above the bound state lies a continuum with energies labeled by the 
wavevector q e (0,7r): E^ — Ej^f + i J — 2<cosg. For every q there are two degenerate wavefunctions ■i/'^ classified 
according to their parity. Odd states are unaffected by the perturbation (5/j_o in Eq- ( ^-21) because the wavefunction 
vanishes on site while even states, which include the bound state, have an on site probability given by: 



l<(o)P = 



J 



32^2 sin^ q 
L [ J2 + I6t2 sin^ q] 



(2.3) 



The spectral function easily follows from the aforementioned properties of the eigenfunctions. In fact, the Lehmann 
decomposition of Aij), uj) gives: 



A(p,w) = <s,~p\cpi\M > \'^5{uj~E,) 



(2.4) 



where the sum is over all the one hole states |s,— p > of momentum —p (modulo tt due to the doubling of the 
cell induced by the antiferromagnetic ordering in \ M > ). These states can be quite generally written as Bloch 
superpositions of eigenstates \s > of the Hamiltonian (2.1): 
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\.s,-p>=^Y.^''"T''\'> (2-5) 



where is the R sites translation operator. This identity readily gives the matrix element appearing in Eq. (2.4) 
in terms of the on site value of the eigenfunction V's(?') corresponding to \s >: < s, — pjcpJA/" >= ^ips{0)- Therefore 
both the energy levels and the matrix elements are independent of p. This is due to the absence of fluctuations in 
the classical Neel state which leads to a local Green function, as already noticed by Brinkman and Rice in the 



J — » limit of this problem. By substituting this result into Eq. (2^) and taking into account the form of the energy 
spectrum we find that, in the thermodynamic limit, the bound state is singled out because of its finite value of the 
on site probability and gives rise to a delta function peak. Instead the other states merge into an incoherent band: 



where lo is now measured from the reference value Ej\/ ~^ \J ^^^d Q{x) is the step function. The incoherent portion 
has several interesting properties: it is a regular, even function of lo which vanishes at band edges oj — ±2t for every 
nonzero value of J. At large J {J > At) it shows a broad maximum at a; = while for J < 4t it has a minimum at 
w = and two symmetrical maxima appear at uj — ±^\/16t'^ — J^. In the J — *■ limit the incoherent part develops 
square root singularities at band edges in agreement with the BR analysis. The shape of the spectral function for two 
representative values of the coupling J is shown in Fig. 1. 

In conclusion: The exact Green function of this problem is purely local and then the hole does not propagate in 
the Ising antiferromagnet. This is conventionally understood on the basis of the "string" defect that the hole creates 
in the antiferromagnetic ordering when it hops [y . However, the quasiparticle weight Z is finite at all non zero values 
of J: The hole behaves as a free particle of infinite mass. This is due to the nature of the excitation spectrum of 
the Ising model which does not allow for gapless modes. In fact we will show that both features of A{p, oj) will be 
strongly modified in more realistic models of hole dynamics. 



III. THE EFFECTIVE HOLE HAMILTONIAN 



In the following we consider the t ~ J Hamiltonian (1.3) defined on L sites with periodic boundary conditions. Our 
task is to derive an effective spin Hamiltonian describing how the hole hopping processes perturb the antiferromagnetic 
background in the particular case of single hole doping. In classical physics this would correspond to a Galileo 
transformation from the laboratory frame to the reference frame locally at rest with respect to the hole. This 
transformation can be easily generalized to quantum mechanics and the derivation can be performed in arbitrary 
dimension. As a first step we notice that the t — J Hamiltonian (1.3) is translationally invariant and then any one 
hole state with definite momentum k and spin | can be written as: 



L-l 
i?,=0 



'ikRrpR 



(3.1) 



where \ao > is a suitable spin state with the spin at the origin R = fixed to | and is the spin translation operator 
defined by the transformation property: 



S 



R+l- 



(3.2) 



where periodic boundary conditions (PBC) over the L sites are assumed in order to define the effect of translation at 
the rightmost site. The latter relation determines the unitary operator Tl only up to an arbitrary phase factor. In 
order to fix the phase of Tl it is enough to specify the action of the operator on a reference ("vacuum" ) state >. 
Here we follow the convention to impose that the ferromagnetic state > is translationally invariant: Ti^\F >= \F >. 
By substituting the representation (3J) of the one hole state into the eigenvalue equation for the t — J Hamiltonian 
we find that > is an exact eigenstate of the t — J Hamiltonian if and only if |cro > is an eigenstate of the following 
effective spin model: 



t [e-'^Ti 



h.c] + J 



Sr ■ Sr+i 



.R=i 



(3.3) 
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where the hole-translation operator Te is defined exactly as (3-2) but with PBC on a squeezed chain of i = L ~ 1 
sites, without the origin i? = where the hole sits. Hk is indeed defined on i sites since the hole at the origin 
is decoupled from the other sites, namely Hk commutes with the spin operator 5*0 at the origin. In the following 



we take the convention to set the spin at the origin with down orientation: Sq 



Notice that the magnetic 



part, proportional to J , represents a Heisenberg model with open boundary condition, since all the magnetic bonds 
connecting the spins with the hole are obviously suppressed. Instead the hole kinetic term K (i.e. the first term in 
Hk) is written in terms of the translation operator which enforces periodic boundary condition on the squeezed 
chain. In this way we effectively traced out, with no approximations, the charge degree of freedom reducing the one 
hole problem to a purely, non translationally invariant, spin model. The effective Hamiltonian Hk explicitly depends 
on the momentum k of the state showing that the distortion of the antiferromagnetic ordering does depend in a non 
trivial (and non local) way on hole motion. 

Also the hole dynamics can be conveniently expressed in terms of the eigenstates of Hk- The spectral function of 
a J, hole is written, in Lehmann representation, as 



A{p,uj) ^ -lTnG{p,uj) 



(3.4) 



where | ^' > is the ground state of the model with no holes and | s > represents a complete set of one hole intermediate 
states. The corresponding energies are respectively Eq and Eg while the momentum space annihilation operator is 
defined by 



E 



(3.5) 



Note that with the adopted definitions, the Heisenberg ground state | ^I* > has total momentum iV7r, where N is 
the number of spins up in the squeezed chain N = Sz + ^/2, and 5*2 is the z— component of the total spin. As a 
consequence, the intermediate states | s > mus t have momentum Nt: — p. By using the general representation (3.1) 
of one hole states with momentum fc, Eq. (|]J) becomes 



A{p, u) 



< cr, |nox|* > \^ d{LO Es + Eo) 



(3.6) 



where now the sum runs over all the eigenstates | (Ts > of H^ with k = Nt: — p. This expression shows that the 
quasiparticle weight of the hole is simply expressed as the modulus square of the overlap between the Heisenberg 
ground state and eigenstates of the effective hole Hamiltonian. 

Equation (3.3), specialized to the case of a hole of momentum p, can be equivalently written in a form which 
makes explicit connection with the problem of a Heisenberg model with a local perturbation. In fact, by adding 
and subtracting the additional bond operator J — J Si ■ Si, the magnetic part of the Hamiltonian H^ can be made 
translationally invariant (on the squeezed chain): 



H„^K + H„-J 



t{-l)^ [e'PTe+ /i.e.] + J 



■i+i 



.4=1 



J Si ■ Si 



(3.7) 



The only extensive, i.e. 0{L), term in the total hole Hamiltonian Hp is H„ which coincides with the usual Heisenberg 
model, while the other two contributions induce 0(1) corrections to the total energy. Therefore, the physics of the 
hulk of the spin system is not affected by the presence of the hole which, at most, acts as a boundary term on the 
Heisenberg antiferromagnet. This can be simply understood in two limiting cases: for a static hole, since t = Q the 
Hamiltonian Hp becomes a Heisenberg model with open boundary conditions, while, for J 0, the last term J is 

irrelevant and the eigenstates of Hp are those of a Heisenberg model with periodic boundary conditions on a £ site 
chain. A close connection between the presence of the hole and a change in boundary conditions on the magnetic 
Hamiltonian thus emerges quite naturally, in every dimensions, within this formalism. 

In one dimension we can proceed further by mapping this spin ^ model into a spinless fermion Hamiltonian via 
standard Jordan- Wigner transformation It is convenient to introduce an additional phase factor in the fermion 
creation operators t/iJj in order to obtain the usual sign of the kinetic energy term: 



S„ 



(-l)'Vtexp 



(3., 
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The magnetic part Her of the Hamiltonian (3.7) becomes: 



H^ = 



I i 

^(^j^.+i + h.c.) + j^(i - 4yj,)il - V'l+iV^.+i) 



(3.9) 



i=l 



where the boundary conditions in the first term are periodic or antiperiodic for even or odd values of L + A^, i.e. 
ipe+i = (— l)^'''^V'i- As usual, N is the number of fermions, which is related to the total magnetization of the original 



spin model by Sz 



N ■ 



Analogously, the bond term J is written as: 



J 



(-1)^+^ + h.c.) + JQ - i^ltPi)Q - ^li^i) 



1 



(3.10) 



In order to express the hole kinetic term K in terms of the spinless fermion operators we have to relate the spin 
translation operator to the usual fermion translation operator Tf which leaves invariant the fermionic vacuum 
state and satisfies TfipiTj = V-'i+i- Keeping track of the phase factors in the definitions we have Ti = {—l)^Tf. 
The fermion translation operator is then conveniently expressed in terms of the Fourier transformed operators ipk = 
£-1/2 ^ ipje^'^^ where the momenta k are quantized according to the choice of boundary condition: 



T, = (~1)^^T/ =exp 



(3.11) 



By use of Eqs. (3. £ ,3.11) we can express the Hamiltonian of a spin down hole of momentum p in the t — J model, Hp 
(3.7), in terms of spinless fermion operators: 



Hp = t exp 



i=l 1=1 



(3.12) 



This concludes our formal manipulations on the original problem. We now have an interacting fermion system which 
can be studied by means of the powerful techniques developed in the framework of one dimensional physics, ranging 
from renormalization group approaches to bosonization methods fl^ . However, before addressing these issues, it is 
instructive to dwell on the similarities between the problem of hole motion and the effect s of local perturbations in 
Luttingcr liquids, which emerge naturally from the structure of the effective Hamiltonian (3.7). 



IV. RELATIONSHIP TO THE IMPURITY PROBLEM 



In this Section we like to investigate the close relationship between the physics of a single hole in the t — J model and 
the impurity problem in Luttingcr liquids. The connection between these two different problems becomes apparent 
when we take advantage of the previously discussed Galileo transformation. The Hamiltonian Hp governing the 
dynamics of the hole represents a Heisenberg model with open boundary conditions plus the hole kinetic contribution 
which involves the translation operator Tj. Furthermore, note that the matrix elements appearing in the spectral 



function (3.6) are related, through Eq. (3.1), to the overlaps between the ground state | G > of the unperturbed 
Hamiltonian {J = K = 0) on L sites and that of the perturbed one (J ^ 0) in the squeezed chain of £ sites. 
The difference in the number of sites of the chain defining the perturbed state can also be thought of as the local 
perturbation induced by the removal of the two bonds connecting the origin in the L site ring. 

Let us focus our attention on the Hamiltonian of a hole of momentum p in the L-sitc t—J model with XY spin 
anisotropy: the t — Jxy chain. In fermion representation Hp is: 



Hr, 



t exp 



j=i 



+1 



h.c. 



N+L 



tpltpl + h.c. 



(4.1) 



where J' — J represents the perturbation induced by the hole. In order to be specific, we consider the model with 
N — 2i/ + 1 fermions (spinous) in a lattice oi L — 1 = £ sites. The kinetic operator K commutes with the magnetic 
term in Ha- and therefore, at J' = 0, the eigenvectors of Hp coincide with those of the XY model on £ sites with PBC 
while the eigenvalues are shifted by 2 1 cos(p + Q) where Q is the momentum of the spin state. The physics of hole 
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motion is contained in the perturbation term proportional to J', therefore it is convenient to study the effects of the 
local perturbation as a function of its strength J'. In the static limit {t = 0) the hole kinetic term K is suppressed and 
we recover the Hamiltonian of a XY model with a weak bond. This problem has been extensively studied in the past 



1^,^: the last term in Eq. (11) is a relevant perturbation and at long wavelengths the system behaves as a chain 
with a missing bond. As a consequence, the overlap ^ between the state with J' = and the state with finite J' goes 
to zero in the thermodynamic limit with a universal exponent X independent of the strength of the perturbation: 

C=<J'\G>(xL-^ (4.2) 

According to conformal field theory, the exact value of the exponent X can be found analytically by taking the 
difference between the size corrections of the open and the periodic chain and dividing by 27tVs, where Vg is the spinon 
Fermi velocity. In the XY model we get X = ^^/4 + 1/16 where /i S (—1/2, 1/2) is the magnetization of the XY 
model along the z direction. 

The recoil of the hole, embodied in the hole kinetic term K, qualitatively changes this picture as we already 
anticipated. The perturbation J' is now marginal: it does not drive the model towards open boundary conditions but 
instead it changes the boundary conditions of the fermionic model in Eq. ( [4.l[ ) introducing phase shifts for right and 
left moving spinous which depend on the strength of the perturbation J' and vanish as J' ^ 0. 

The effect of the kinetic term can be understood by use of first order perturbation theory in the parameter J' in 



Eq. (4.1). For J' = the ground state |G > is non degenerate and can be represented by a Slater determinant of 



plane waves with PBC. At first order, the ground state becomes: 

1^' >- |G > +^ E + ^-'') ^'''^'^^^^"'^ \k, 1 > (4.3) 

^'^ k^q ^^^1 

where the sum runs over the momenta of particlc-hole excitations of the unperturbed system and e^^q are the corre- 
sponding excitation energies: 

Cfc.g — J [cos q — cos k] + 2t [cos(p + fc — q) — cosp] (4.4) 



The matrix element in Eq. (4.2) < A;, gj-^I^^AqlG > is unity provided q belongs to the Fermi sea while k lies outside 



the Fermi surface. Perturbation theory fails in the thermodynamic limit due to low energy excitations which may be 
of "forward" (i.e. q ~ ±.kp and k ~ ±fci?) or "backward" type (i.e. q ~ ±.kp and k ^ T^f)- However, for any t 
the backward scattering is cut off by the recoil and does not introduce singularities in perturbation theory, as can 



be easily c heck ed by use of Eq. (4.4). Instead, forward scattering is always singular leading to the vanishing of the 
overlap C (4.2) in the thermodynamic limit. 

In order to better understand the change in the state induced by the boundary term J', let us consider a different 
problem, namely, the way the ground state changes due to a weak modification of the boundary conditions. Let 
us take as unperturbed state a Slater determinant of plane waves 4>q{r) — e^'^^l^^/^ with momentum quantization 
appropriate to PBC: £qj — 27rj. A change in boundary conditions involves the introduction of certain momentum 
dependent phase shifts dk- Therefore, we take as perturbed state another Slater determinant of plane waves ^k{f) 
with generalized momentum quantization rule: tkj = 2-K{j + 5k)- Here and in the following, the phase shifts will be 
measured in units of 2-k. For weak perturbations, the phase shifts will be small and we can keep only 0{5k) terms in 
the expansion of the state in the unperturbed basis. As usual, to linear order, only a single particle-hole excitation 
|fc, g > is allowed and the weight of the corresponding contribution turns out to be: 



(4.5) 



for k ^ q. Again, we see that a change in boundary conditions leads to singularities only for forward scattering as 
in the hole problem when the effect of the recoil is taken into account. T he com parison can be made quantitative by 
linearizing the momentum dependence about the Fermi points in Eqs. (4.5,1^) and matching t he t wo expressions. 



The effective phase shifts induced by the perturbation J' at long wavelength in the Hamiltonian (4T) are then given 
to O(J') by 

2-,5^ = /^^"^^^ , (4.6) 
2t sinp ^ J sin kp 

at the Fermi points. The two signs refer to the right and left moving spinous which in general have different phase 
shifts. Notice that the resulting phase shifts depend, in general, on the density as well as on the total momentum 
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p of the hole. Following the well known analysis of Anderson's orthogonality catastrophe ||21|], we then find that 
the overlap between the states before and after the perturbation should vanish in the thermodynamic limit with an 
exponent X = (S"^ + (5^)/2. Few details are reported in appendix (^. 

In conclusion, wc can interpret the effects of the presence of a hole in an antiferromagnetic background as a change of 
boundary conditions in the corresponding spin problem. If the recoil of the hole is neglected, (i.e. t — 0) both forward 
and backward scattering between the holon and the spinous are relevant. In this case, the hole effectively breaks the 
spin ring and at low energy the model becomes equivalent to an open spin chain (in the hole reference frame). Instead, 
for t ^ the backward scattering channel is cut off and the perturbation induced by the hole becomes marginal. The 
presence of the holon introduces, via the forward scattering channel, phase shifts in the boundary conditions of the 
spin chain. Spinous can then propagate through the site where the hole sits. 

These results have been obtained analytically for a XY model but, in view of the universality of the Luttinger 
Liquid description of ID correlated systems, we expect that our picture remains valid for general XXZ spin chains, 
including the isotropic Hcisenberg point. 



V. FIELD THEORETICAL ANALYSIS 

In this Section we address the problem defined by the Hamiltonian Hp ( 3.12| ) in fermionic representation by use 



of the bosonization method. The bulk term H^r (3.9) just represents a Heisenberg model with periodic boundary 
conditions which has been extensively studied in the past The renormalization group approach applied to 

such a problem shows that both backward scattering and Umklapp scattering terms are marginally irrelevant: at 
long wavelength the model is characterized only by the interactions of the forward type with non universal coupling 
constants renormalized by the RG flow. This fixed point Hamiltonian represents a Luttinger model which can be 
exactly mapped into a free bosonic system following the procedure of Mattis and Lieb , the bosons representing 
density fluctuations. As suggested by the analysis of Section the additional terms present in Hp slightly modify this 
picture, appropriate for the undoped antiferromagnet. The combined effect of K and J on the long wavelength physics 
is equivalent to a change in the boundary conditions of the Luttinger model from periodic to skew: Although the 



original problem (3/7)is not translationally invariant due to the bond term J, the low energy fixed point Hamiltonian 
develops effective boundary conditions which depend on the non universal parameters characterizing the microscopic 
model. In particular we expect that the boundary conditions will continuously evolve from periodic at small J/t to 
generally skew at finite J/t with the exception of the zero magnetic field case where the boundary conditions at the 
fixed point remain periodic for any t ^ due to the presence of the additional SU (2) spin symmetry, as will be shown 
later. The only singular point is at t = where the hole kinetic contribution vanishes and Hp describes a Heisenberg 
model with open boundary conditions. In the following, we will exploit this picture extracting quantitative predictions 
which will be later compared with the exact solution in specific models. 

On this basis we are led to consider a low energy problem defined by a long wavelength effective Hamiltonian Hp 
sum of two commuting terms: a Luttinger model with skew boundary conditions plus the hole kinetic term 

K = ehikh) (5.1) 
kh^P + Q (5.2) 

Here we have defined as the holon momentum and Q is the spinon momentum which obeys the quantization rule 



appropriate for the skew boundary conditions of H^. Eq. (5.2) therefore represents momentum conservation of the 
charge and spin excitations. Generally, we expect that the renormalization group flow modifies the effective hole band 
in a non universal way: 2tcos{kh) th{kh)- 

Clearly, such a form of the long wavelength Hamiltonian may describe only the low energy part of the excitation 



spectrum of the effective hole Hamiltonian (3.12). Therefore, an elementary spin excitation which changes the total 
spinon momentum Q must be accompanied by a corresponding variation of the total momentum p so that the holon 
momentum kh is not modified. Otherwise the hole kinetic term K would contribute a finite amount to the excitation 
energy and a low energy excitation in the spin channel would not correspond to a low lying excitation of the total 
Hamiltonian. In other words, the low energy spin excitations described by the effective spinon Hamiltonian with skew 
boundary conditions represent the physical excitation spectrum at fixed holon momentum k^. 

The overall picture of the low energy physics of the one hole problem is confirmed by the exact solution of Bethe 
Ansatz models which also allows for a quantitative analysis of their long wavelength properties as will be shown later. 

As a first step, we now generalize the bosonization procedure in Luttinger models with skew boundary conditions 



(Sec. VA), then we analyze in some detail the form of the hole kinetic operator present in the effective Hamiltonian, 



the possible boundary conditions which may occur in microscopic models and the corresponding finite size corrections 
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to the ground state energy of this effective model (Sec. VB). Finally we discuss some i mplic ation of the previously 
obtained results including the expected asymptotic form of the hole Green function (Sec. V_C). 



A. The long wavelength Hamiltonian: Bosonization 



At low energy, the important degrees of freedom for a many fermion system are those close to the Fermi momenta 
kp for the left and right movers, i.e. k ~ +kp in the right (+) branch and k ~ —kp in the left (— ) branch. As usual, 
the two branches are extended to infinity, within the assumption that the low energy physics is not affected by this 
approximation ||l3|| . This extension allows to define two fields ip+i^) ^-nd representing, in the continuum limit, 

the annihilation operators for spinless fermions on the right and left branch. The continuum limit of the original 
fermionic field ip''{x) is then given by the linear combination: 

^-^(2;) cx e^'^^il^lix) + e^^'^^^i^lix) (5.3) 

where we have kept the distinction between right and left Fermi momenta kp which may in principle differ. The 
origin of the fermionic field x = has been chosen to match with the first site z = 1 of the squeezed chain of length 
£. As stated before, the fields ip±{x) obey skew boundary conditions: 

i'±{x + e^'"^±i;±{x) (5.4) 

with arbitrary phase shifts S± . 

The Luttinger model is defined by right and left moving fermions with kinetic term and interactions of the for- 
ward type. This allows to express the effective low energy theory in terms of two bosonic fields defined as bilin- 
ear combinations of the fermionic operators. More precisely, following Mattis and Lieb the Fourier transform 
-^±(9) — J dxe~^'^'^ N±(x) of the operators N±{x) — i^\.{x)t}:±{x) + const, satisfies non trivial commutation rules 

[N±{q),N±{-q)]^T^ (5.5) 
The above equations suggest the definition of the boson field $(a;) and its conjugate momentum n(a;) via: 

N+{x) = i^l{x)i;+{x)- < ijl{x)^+{x) > = {Il{x) + d,^{x)) 



N-{x) ^ipl{x)i:-{x)- <^pl{x)i:-{x) > = = (n(x) - a^$(x)) (5.6) 



where the average < >, taken on the reference ground state of the Luttinger model, is introduced in order to regularize 
the divergences. These operators obey the canonical commutation relations: 

[<i>{x),Il{x')]^iS{x-x') (5.7) 

Notice that the densities N±{x) satisfy periodic boundary conditions independently of the phase shifts 6±, being 
bilinear combinations of fermionic operators. As a consequence the fields dx^{x) and Tl{x) obey periodic boundary 
conditions and the Luttinger Hamiltonian, which can be entirely expressed in terms of the bosonic fields ($(a;),n(a;)), 
has exactly the same gaussian structure as for a pure Heisenberg model: 



-f^a = y J dx 




const (5.8) 



where Vs and are the renormalized Fermi velocity and the dimensionless interaction parameter which characterize 
the long wavelength behavior of the Heisenberg model . 

This simple quadratic Hamiltonian can be diagonalized by introducing normal modes and this procedure leads 
to the familiar Luttinger Liquid energy spectrum of the one dimensional Heisenberg model. However, here we are 
interested in the change in the finite size corrections to the energy spectrum induced by the presence of the hole. 
Therefore, the normal modes must be carefully defined in the £ site chain. As a first step note that two quantum 
numbers labeling the eigenstates can be defined: The two operators 
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A$ = [$(^) - $(0)] ^ Vt^ W+ix) + N^{x)] dx 



(5.9) 



commute each other and also commute with Ha-, being related to the total number of fermions on the two b ran ches, 
which are conserved quantities in the Luttinger model. Therefore the diagonalization of the quadratic form ( |5.8| ) can 
be performed in each sector defined by the pair of quantum numbers (A$, 11*). In particular, the choice A$ = 11* = 
identifies the reference state introduced in Eq. (|5.6| ) . The quantization rules directly follow from the definition ( ^.9| ) of 
these operators. At fixed total number of fermions (i.e. at fixed magnetization in the original model) A<i> is uniquely 
determined. A low energy excitation in this model corresponds to moving a fermion from the left to the right branch. 
This changes the value of H* by an even multiple of ^/n. On the other hand, by varying the total magnetization 
by integers, A<I> changes by integer multiples of ^/tt while the model preserves a non degenerate (degenerate) ground 
state for even (odd) £ due to a change in the boundary conditions associated with the Jordan- Wigner transformation. 
Therefore we conclude that the possible values of A<i> and H* are 



A$ 

n* 



2\/Tr{m 



7) 



(5.10) 



Here m and n are integers while 7 and S are two non universal real quantities (defined mod(l)) which characterize 
the ground state of the model. Their precise value is determined by the renormalization group flow which connects 
the microscopic model to the Luttinger Hamiltonian. The bosonization procedure alone does not fix these quantities 
uniquely, except when additional symmetries are present in the microscopic model. 

Now we are ready to introduce the normal modes which diagonalize the quadratic form (5.8). It is convenient to 
define the periodic field 



By substituting into (5.8) we get: 



"ifix) = $(x) - -A$ 



(5.11) 



Ha = ^ I dx 



21K, 



-(A$)2 



(5.12) 



Next we define canonical creation and annihilation boson operators: 



+ a-fc) 



dxY\{x)e 



ikx 



2Ka 



(5.13) 



where the values of fc 7^ correspond to periodic boundary conditions and are quantized in units of Finally we 
get the low energy Hamiltonian in diagonal form: 



21 \ Ka 



(5.14) 



Here Eq is the reference energy of the non degenerate ground state of the model with 7 = (5 = which is known, from 
conformal field theory H, to depend on the lattice size i as: 



£■0 = eo^ - Vs 



(5.15) 



if Co is the Heisenberg ground state energy per sit e. N otice that whenever 5 or 7 acquire non zero values only the 1 / i 
finite size corrections to the ground state energy ( [3.12| ) are modified. The quantum numbers 7 and 5 then correspond 
to energy contributions induced by a change in the boundary conditions of the model. This remark will be made more 
precise in the following. 
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B. Spinon momentum and skew boundary conditions 



In the previous Section we investigated the bosonization of the bulk part of the spin Hamiltonian in some detail 
At t he Luttinger fixed point the Hamiltonian H^r can be expressed i n te rms of bosonic operators in the standard form 
(5.8) independently of the presence of skew boundary conditions (5^)- However, the energy spectrum depends on 
two quantum numbers (A<I>,H*) which should be somehow related to the particular choice of boundary conditions. 
In order to make this relationship more transparent, let us analyze the bosonization form of the momentum operator 
(i.e. of the translation operator) appropriate to the chosen boundary conditions. 

The fermionic translation operator is defined by Tfi/j{x)Tj = ip{x + 1). This unitary transformation can be split 
conveniently into two steps: first a conventional translation of left and right movers '>jj±{x) ~^ '>jj±{x + 1) (5.17) and 
the n a rescaling of the fields: ?/'+ ~^ e"*'^"^ 



^■0+ and "0- e 
in (5.3). The correct form of Tf is then Tf — T1T2 where 



Fip- ( 5.16 ), in order to take into account the phase factors 



To 



exp 



exp 




k+ipl_{x)ip+{x) - kpi)'^_{x)i^-{x) 



ipUx)dxi^+{x) + ■4;^_{x)dxik-{x) 



(5.16) 



(5.17) 



There is no overall constant in the definition of the translation operator since the e xpre ssion (5.17) acts as the identity 
on the fermionic vacuum state. Note that Ti has been already bosonized via the (5.6). The other factor T2 is just a 
conventional translation of the fields 4'±{x) f-nd ^^cts on the bilinear forms N±(x) as the bosonic translation operator: 
H(a;) ^ T2H(a;)r2"^ = U{x + 1) and $(a;) ^{x + 1). Therefore 



Ta = exp I ^ y" dx [U{x) d,<i> + d^<i> n{x)] | 



(5.18) 



From the previous analysis is clear how to represen t th e translation operator Tf within bosonization. By use of the 
explicit expressions of the terms in Ti , T2 and Eq. (|5.6|) we get 



T 



f 



P = Qo 



20F 



-n* 



2y/^ 2 



dx [n{x)d^<i> + d^<S>U{x)] 



Qo 



-{n + -/) + ^J dx [Uix) + U{x)] 



(5.19) 



where in the last line we took advantage of the canonical transformation (5.11) and we expressed H* and A<I> by use 
of Eqs. ( ^.10 ). The additive constant Qo represents a finite contribution which depends on the ban dwidth cut -off o f 
the Luttinger model but not on the chain length £ and is obtained by the substitution of Eq. (5.6) into Eq. (5.16). 
Its value is determined by requiring that in the ground state (n = m = 0) the spinon momentum, in the £ ^ 00 limit, 
tends to kp, i.e. t o th e momentum of the translationally invariant case without perturbation induced by the hole 
( J = = in Eq. |3.7| ). This gives Qo = kpi^ — 2S) for odd £. The total momentum operator P commutes with the 
bulk Hamiltonian and then each spinon state is characterized by a spinon momentum Q which is related to the 
total momentum p and to the holon momentum kh by Eq. (5.2). 

Through the bosonized form of the fermion translation operator it is also possible to relate the phase shifts S± 
which define the boundary condition of the microscopic model to the effective parameters (7,<5) entering the low 
energy Hamiltonian. The elementary excitation which changes the value of H* corresponds to moving a fermion from 
the left to the right branch. In fermionic representation, this excitation carries momentum kp+kp + ^{S+ — 5-). 



Instead, in bosonic representation, this corresponds to n 
momentum due to such an excitation we have to identify 



and 



1. In order to match the change in total 



7 = (5+ — (5_ mod(l) 



(5.20) 
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Analogously, a low energy and low momentum excitation which changes the total spin corresponds to adding one 
fermion to the right branch and one fermion to the left branch. The change in momentum is then kp — kp + '^{S^+6-) 
while the quantum numbers are n — 2 and m — 0. Matching the two expressions gives: 

6 = ^^-^ mod(l) (5.21) 

This concludes the bosonization of the long wavelength effective Hamiltonian, which includes the bulk contribution 
Her and the hole kinetic term K which is written in terms of Tf. Collecting the various terms together our final 
expression for the low energy spectrum of one hole at momentum p is: 



Ep 


— E^ + Eh 






(5.22) 


Ea 




I 




(5.23) 


Eh 


= Qiikh) 






(5.24) 


kh 


^p + Q 






(5.25) 


Q 


= kp (2m + 1) 


27r r 
+ -[ 


(tj+ + n + 7)(r7i + (5) + -^(" + 7) 


(5.26) 






(mod 1). In the thermodynamic limit, the quantization rule of the spinon 



momentum reduces to the non interacting result which, for odd £, reads 

Q,„ = fcF(2m + l) (5.27) 

but the finite size corrections of energy and momentum explicitly depend on the phase shifts. We stress that both 
the parameters entering the effective, long wavelength Hamiltonian (ws, Ka, kp) and the phase shifts {S±) cannot 
be trivially related to the bare lattice Hamiltonian because the RG flow renormalizes all the couplings not protected 
by conservation laws. Their value can however be uniquely determined in Bethc Ansatz soluble model by matching 
the form of the finite size corrections of the one hole energy and momentum. This program will be pursued in the 
following Sections. More information can be gained at zero magnetization. In this case, the effective spin model ( |3.3| ) 
has the additional SU(2) spin rotational symmetry which limits the possible boundary conditions of the Luttinger 
model. In fact, in spin isotropic models, the allowed boundary conditions are either open or periodic: non trivial phase 
shifts are not compatible with the requirement of spin isotropy. Therefore, according to our basic assumption the hole 
kinetic term stabilizes the periodic boundary conditions in the effective spin Hamiltonian. This observation leads to 
a unique determination of the phase shifts S± at zero magnetization which only depend on the parity of the number 
of sites i. We first note that the z component of the total spin is simply related to the total number of fermions 



and then the spin excitations are labeled by the quantum number n in (5.1C). By definition n = corresponds to 
the ground state. Moreover, in zero magnetic field, states with opposite values of are degenerate. For even £, the 
groun d state is a singlet and then excited states labeled by n and —n are degenerate: This implies 7 = through Eq. 



( 5.23 ). For odd £ the ground state itself is a spin doublet which gives 7 — 1/2 (mod 1) by the same argument. For 



periodic boundary conditions also the total momentum is a g ood quantum number. At fixed magnetization, excited 



states of definite momentum are labeled by m in Eq. (5.10), the ground state corresponding to the choice m — 0. 



For e ven £ the ground state is unique and excited states with opposite momentum are degenerate giving 6 — hy 



a 



(5.23) while, for odd £ the ground state has finite momentum which implies a twofold degeneracy due to parity. Such 
degeneracy is compatible with the form of the energy spectrum ( [3.23| ) only for S = 1/2. 
As an example of the previous analysis, we plot in Fig. 2 the overlap square Z between the Heisenberg ground state 
on a £ = L — 1 site ring and the ground state of a hole of momentum p = 7r/2 in the L site t — J model at J = At and 
vanishing magnetization. The results have been obtained by Lanczos diagonalizations in chains with even L < 26. 
According to our analysis we expect that at long wavelengths the single hole problem in an even chain is described 
by an effective Heisenberg Hamiltonian with (5 = 7 = 1/2 which correspond to periodic boundary conditions on i — 1 
sites, leading to a finite overlap in the thermodynamic limit. Actually, the numerical results provide quite a strong 
evidence in favor of this picture showing an overlap Z which increases with the size of the system. 

In conclusion we see that spin isotropy determines the quantization constant 7 while the occurrence of periodic 
boundary conditions in the effective long wavelength Hamiltonian fixes the value of S. At finite magnetization, spin 
isotropy is broken and we expect the occurrence of generic, momentum dependent values for 7 and 6. These phase 
shifts will play an important role in determining the singularities of the Green function. This subject will be discussed 
in the following Section. 

For completeness let us briefly extend the previous discussion of the energy spectrum of the model to the case of 
open boundary conditions. According to our assumption this case is relevant only when the hole effective hopping 
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amplitude t vanishes. The general form of the bulk Hamiltonian 7?^ (5.8) does not depend on the choice of boun dary 
conditions and is therefore unaltered. However, after having performed the canonical transformation ( |S.ll ), the 
normal modes are now defined by use of standing waves rather than the previously introduced plane waves (5.13): 



2 



sin(fc„a;)n(a;) = ^-(^l^ - ^kj 



where now fc„ 



(2rt+l)7r 
2£ 



with n > 0. The bulk Hamiltonian for open boundary conditions then reads: 



Eq 



n>0 



■(n + 7)2 



(5.28) 



(5.29) 



The final expression for the energy spectrum of the hole problem at t = 0, i.e. when open boundary conditions apply, 
is then 



E — En 



2nv^ 1 



where now the size scaling of the reference energy is p9| 



■(n + 7)2 



(5.30) 



(5.31) 



Again, in the spin isotropic case the values of the phase shifts are constrained by the SU{2) symmetry which gives 
7 = for even I and 7 = 1/2 for odd t. 

Now we conclude this Section by expressing the expected finite size corrections to the one hole energy which emerge 
from our picture. On the basis of the discussion at the beginning of this Section, it is particularly convenient to work 
at fixed holon momentum kh extracting the size corrections of the low energy excitation spectrum obtained by varying 
the spinon momentum and the magnetization of the model. These results will be later compared with the exact form 
of the energy in two Bethe ansatz soluble models. The size dependence of the hole kinetic contribution (5.24) can be 
obtained by direct substitution of the explicit form of the spinon momentum ( 5.26| ) giving a term proportional to the 
charge velocity 



Vc(kh) 



dkh 



(5.32) 



Here kh is related to Q,„ and to the total momentum p by the conservation law ( ^.25 ). We also allow for a size 
dependence of the total hole momentum p: 



P^Po + 



2'Ka 



(5.33) 



valid up to o{l/i) terms. The constant a depends on the adopted sequence of lattice sizes. By the conservation law 
kh = P + Q (Eq. ^.25 ) the holon momentum acquires the finite size corrections of p (Eq. 5.33) and Q (Eq. p3.26 ). 
Therefore, the holon kinetic term (5.24) will contribute to the 0{l/t) size corrections of the ground state energy as 



2'KVr 



a + ((5 + m) (7 + n) + w^. (J + m) + -w_ (7 + n) 



(5.34) 



On the other hand, the spinon term, being described by the conformal fi eld t heory which characterizes the Luttinger 
liquids, gives rises to the finite size corrections already obtained in Eq. ( 5.23| ). 

In conclusion, at fixed holon momentum fc/j, the finite size corrections of the energy Ep of the single hole problem 
give rise to a tower of states which depends on two quantum numbers (n, to) and have total momentum p = kh — Qm 
(5.27). The size corrections are uniquely determined by three bulk properties of the system (wc, Ws, K^) and four 
additional constants which determine the boundary conditions of the effective spin Hamiltonian (w±, 7, (5): 



AEr, 



2TTVr 



2'KVs 



a + (^ + to) (7 + n) + w+ ((5 + to) + -tj- (7 + n) 



(5.35) 
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Through this equation we can in principle evaluate the elusive spinon phase shifts S± by computing the easily accessible 
finite size corrections of the one hole energy. This equation has been obtained in the framework of the t — J model 
where no double occupancy is allowed but we will show that the same structure persists also in the one dimensional 
Hubbard model at finite U. Therefore we believe that this form of the finite size scaling is a general feature of one 
dimensional correlated models. 



C. Orthogonality catastrophe and the hole Green function 



In this Section we like to relate the previously introduced phase shifts with the behavior of physical quantities and 
specific dynamical correlation functions of the one hole problem. In particular, we will address first the evaluation 
of the hole quasiparticle weight at holon momentum kh and then the calculation of the asymptotic behavior of the 
Green function. 

The quasiparticle weight is defined as the square of the modulus of the matrix element: 



,1 



* > 



(5.36) 



where \kfi > is the exact one hole ground state at holon momentum /c^, | 5* > is the ground state of the Heisenberg 
model and we have chosen the convention of creating a hole with | spin projection. Momentum conservation implies 
that in the ground state fc/j = p± kp, where kp is the spinon Fermi momentum. By use of the Galileo transformation 
on the hole problem and a further Jordan- Wigner transformation on the "up" spins, the problem is reduced to the 
evaluation of the overlap between the two fermionic states corresponding to the Heisenberg ground state | ^' > on a 
L site ring and the ground state | kh > of the effective spin Hamiltonian defined on the lattice oi £ = L — 1 sites and 
the same number of up spins. In order to study the behavior of the quasiparticle weight for L — > cx), we can limit 
our attention to the long wavelength form of the effective spin Ham iltonian. The problem is therefore to compute 

with different boundary conditions 



the overlap of two eigenstates of the same bosonic Hamiltonian (5 



reference state | vj/ > has standard periodic boundary conditions (i.e. H* 



while the 

— A<i> — 0), the one hole ground state | kh > 



is defined by non universal, momentum dependent, values of the phase shifts: 

A$ = y/n^ = y/n{6+ - (5_) 



(5.37) 



In the continuum limit this overlap is strictly zero because the two states are eigenstates of the bosonic Hamiltonian 
with different eigenvalues for A<i> and H*. This orthogonality is however an artifact of the Luttinger extension to a 
system containing an infinite number of particles and it is clear that in a finite system the overlap between two states 
with different boundary conditions will be in general finite. The solution is to "regularize" the Luttinger branches, 
in order to be consistent with a tight binding model where the same phase shifts 6+ and S- at the Fermi energy are 
induced by a local potential which gives no contribution to the backward scattering. In a free Fermi gas, = 1 
and the two states can be written as Slater determinants of plane waves with suitable phase shifts. In this case, the 
overlap can be easily calculated giving an asymptotic power law behavior which only depends on the value of the 
phase shifts at the Fermi energy of the right and left branches: 



C = < kh\^ > oc exp 
1 



1 



is'. 



(5.38) 
(5.39) 



where the last equality in Eq. ( 5.39| ) follows from Eq.(5.37). The formal calculation of the overlap is contained in 
appendix This formula generalizes the exact result valid for 6+ = 6- [pl] , p3| and agrees with the prediction of 
conformal field theory relating t he fin ite size corrections of the energy to the exponent of the boundary operator s [[l9[ . 
In fact, the explicit expression (5.39) coincides with the term proportional to lnvsj (■ in the final formula Eq. ( |5!35| ) 
for the finite size corrections of the one hole energy. 

This formalism can be easily generalized to the case K^j ^ 1: In fact, the overlap does not change upon unitary 
transformations and it is known that the scaling of the bosonic fields 



H'(a;) = V^n(a;) 
$'(.t) = $(a;) 



(5.40) 
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maps the interacting problem with ^ 1 to the one with = 1. By the same transformation, the boundary 
condition for the one hole state are modified due to Eqs. (5.9, 5.10 ): 7' = ^l\JK„^ 6' = Sy/Ka, while for the reference 
state the boundary condition arc unchanged. The effective phase shifts of the parent non interacting case then follow 
immediately: 




leading to the general expression: 



Xn = 



(5.41) 



(5.42) 



again consistent with the finite size corrections to the energy ( ^.35| ). 

The relation between the finite size corrections to the energy and the exponent Xq of the orthogonality catastrophe 
is a general property of all conformal field theories [|l9| and holds also in other cases. For instance, the exponent 
assoc iated to the open boundary fixed point can be determined in terms of the finite size corrections obtained in Eq. 
(B: 



Xa = 



1 

16 



(5.43) 



where use has been made of the known additional contributions to the ground state energy with periodic (5.15) and 
open (5.31) boundary conditions on a L and £ site chain respectively p9| . This expression, which applies in the limit 
of a static hole, yields the exact exponent Xq = 3/16 in the isotropic case, when = 1/2 and 7 = 1/2. This has 
been checked numerically in Fig. 3 by Lanczos diagonalization of the hole problem in the limit of vanishing hopping 
amplitude t. 

Now we can proceed to the study of the asymptotic behavior of the retarded single hole Green function: 



^-i{H-iri)T^ 



(5.44) 



valid for r > with 77 = 0"*" as convergence factor. As usual, by performing the Galileo transformation and then a 
Jordan- Wigner transformation we can map the problem of the evaluation of the Green function to the calculation of 
the purely fermionic matrix element: 



Gx(p,r)=^ < *|e^*(^''~'")"(l-no)|* > 



(5.45) 
(5.46) 



where the effective Hamiltonian Hp is defined in Eq. ( 3.12| ) and the fermion projection operators ng = ipQipo , 
(1 — no) = ipo4'o force the spin at the origin of the chain to point upwards or downwards respectively. In the following, 
we will first carry out the calculation of the asymptotic behavior of the trace of the Green function matrix: 



G{p,t) = Gi{p,t) + G^{p,t) ^ I < 



-i(Hp-ivi)T 



* > 



(5.47) 



which is expected to show all the singularities present in the two separate spin projections. 

The long wavelength form of the Hamiltonian Hp has been already studied in the previous Section. Hp depends on 
the total momentum p via a function eh{kh) which represents the holon dispersion: 



Hp^eh{p + P)+H„ 



(5.48) 



Hn is a Luttinger liquid Hamiltonian characterized by suitable boundary conditions and commuting with the spinon 
momentum operator P defined on the chain of L — 1 sites. Note that the Heisenberg state | ^ > instead refers to the 
L-site ring with periodic boundary conditions and therefore is not an exact eigenfunction of the spinon momentum 
operator P. It is then convenient to express the state | vjf > as the sum of all its projections into subspaces of definite 
momentum P ~ Q and substitute this representation into the form (5.47): 
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G(p, t) = jY1 ^(Q' r)e"*[^''(P+'3)-"'l " (5.49) 
Q 

Z(Q,t) = / e-''5^Z(x,T)dx 
Jo 

Z(a;, r) = < * | e'^^^^""''^ \ * > (5.50) 

This is the general form of the one hole Green function which shows the effects of spin-charge decoupling on the 
dynamics of the hole. The spinon function Z{Q, t) provides a generalization of the quantity Z{Q) which characterizes 
the form of the single hole Green function in the J ^ limit jllU23]. Here, however, non trivial dynamics of the spins 
are induced by H„ as long as Vs 7^ 0. 

Let us first discuss the calculation of the matrix element Z{x,t) in the non interacting limit. In a free Fermi gas 
with arbitrary boundary conditions, the long wavelength Hamiltonian reduces to the sum of the kinetic terms for the 
two branches of right and left moving fermions. In the continuum limit the Hamiltonian and the momentum operator 
are given by: 

H ^Vs{P+ - P-) + const 
P = P++P- 

P±^±kF dxipl{x)i^±{x) + i / dxiljl{x)da:ip±{x) (5.51) 
Jo Jo 

This particular form of the Hamiltonian shows a close relationship between energy and momentum operators which is 
clearly valid only in the non interacting limit. However, the more general interacting Luttinger liquid can be mapped 



to the free Fermi gas by the previously defined canonical transformation (5.40) leading to the conservation of the 
number of right and left moving fermions at long waveleng th. A s a consequence, all eigenstates would factorize in the 
product of two states, one for each branch and, via Eq. ( ^.50 ) also the function Z(x,t) would split in the product 



of two terms {x =F Vst) defined in each branch. However, in the microscopic model, defined on a lattice, higher 
order terms allow the excitation of a fermion from one branch to another (even for = 1). Therefore, the function 
Z{x, t) will contain contributions corresponding to all these excitations which we label by the number m (—m) of 
extra fermions on the right (left) branch: 

00 

Z{x,t)^ e^(«"--^™^)z|'+'")(a;-t;,T)zi''-'-'"^(x + i>,T) 

m — — 00 

00 

= X! e''^'^"^"-^'"^) < (.iiP+-Q+i"^))i^-^sT) ^i{P--Q-{-m,)){x+v,T) ^_ {5.52) 



where the average labeled by + (— ) is taken on the ground state of the right (left) branch in the undoped system 
and the intermediate states are constrained to have m (—m) additional particles on the right (left) branch. Qm = 
Q +{m) +Q -{—m) and Em = Vs{Q+(m) — Q-{—'m)) represent the reference momentum and energy of the intermediate 
states with m-particle excitations which are explicitly given in terms of the Fermi momentum of the spinous by 



Q+(m) = kp{m + 1) and Q^{—m) — —kp{—m). The functions z£**'^™''(x, r) introduced in Eq. ( ^.52 ) are well 



defin ed also in the long wavel ength li mit w here the spinon Hamiltonian H„ can be written in the bosonized form 
(5.8). As discussed in Sections V A and VB the quantum number m which characterizes the inter-branch excitations 
appears in the low energy Hamiltonian only throug h the quantity H* and can be absorbed in the definition of the 
phase shifts: Sj^ — 5+ + 5- = 6- + m (see Eqs. 5.10, 5.20, 5.21). The calculation of Z{x,t) then reduces to the 
evaluation of the contribution appropriate for each branch with arbitrary phase shift S±. The tec hnica l details are 
discussed in appendix & here we just report the asymptotic behavior in the thermodynamic limit (B.22): 



Zl{x) cx {Tix) 



-5^ 



(5.53) 



By substituting this result into Eq. ( ^.50 ), the spectral weight A{p, uj) ~ ^IraGip, oj) is written in terms of the (Q, uj) 
Fourier transform of Z: 



dQ 

271 



ZiQ,LO~eh{p + Q)) 



(5.54) 



16 



oo /'OO 



Z{Q,uj)= / dTdxe^'^''^-'^^'>Z{x,T) 

J —OO J — oo 

= ^ C„i e [tj + Vs{Q - Q„i)] 9 [w - Vs{Q - Q„i)] 

m 



(5.55) 



In the Fourier transform of Z{x,t) only the singular contributions have been included. Q{x) represents the step 
function whose presence is a direct consequence of the Fourier positivity (negativity) of the function (Z-), as 
discussed in appendix This property has a simple physical meaning: The particle hole excitations within the 
right branch of a Luttingcr model can only increase the total momentum with respect to the ground state, while the 
excitations on the left branch can only decrease the total momentum of the system. Therefore the spinon spectral 
function at Q > has contributions coming only from excitations in the right branch and vice versa. The constants 
Cm, which only depend on the phase shifts (5+ and S- in a sym metri c fashion, can be explicitly calculated in the 
free Fermi gas with skew boundary conditions as shown in Eq. ( B.24 ) of appendix However, we expect that in 
the interacting models these coefficients will be renormalized in a non universal way which depends on the physical 
cut-off present in the mic roscopic mo del. 

By inserting Eq. {j.bl) into Eq. (5.54), the contribution coming from each excited state m gives a divergence in 
the spectral function as long as 5^ + S±< 1. The singularities are located along the lines 



Alu„i = uj - eh{p + Qm) = 
in the (p, uj) plane and show the asymptotic behavior: 



(5.56) 



(5.57) 



where Xm = 5(^+ + ^?.). In general, we expect that the phase shifts 5± which characterize the singularities of 
the spectral function depend on the holon momentum kh- Therefore, A{p,Lo) will show singularities along the lines 
UJ — eh{p + Qm) with momentum dependent exponents X^- 

Another unexpected prediction of the present formalism concerns the behavior of the spectral function in a neigh- 
borhood of the singularity: By a direct evaluation of the prefactor we find different results according to whether 
the charge velocity is larger or smaller than the spin velocity. In fact, if |wc| < Ws, there are divergences only for 
^ ^ ^h{p + Qrn) while the spectral weight vanishes on the other side of the singularity line AtOm = as there are no 
states with lower energy contributing to the spectral weight: 



A{p, uj) cx 



\AuJr 



\2X„ 



(Vs + Vcf+{Vs - VcY- 



(5.58) 



Instead if jwd > Ws the spectral function diverges on both sides of the singularity with the same exponent 2X„ 
but different prefactor: 



A{p,uj) cx 



lAuJr 



\2X„ 



(5.59) 



wher e the upper sign and the lower one refer to zLvcAoj > 0, respectively. The amplitude ratio can be evaluated from 
Eq. (|^: 



A{-Auj) 



sin(7r5^) 



(5.60) 



This simple expression, which shows the asymmetry of the spectral weight above and below the singularity, is strictly 
valid only for the Luttinger model, but is also expected to be qualitatively correct in more general cases. In fact this 
feature originates from the different energy spectrum of the excitations on the left and right branches, which is a 
common property of one dimensional correlated systems. Notice also that the prefactors in ( [5.58 , 5.59) are strongly 
enhanced close to the instability which sets in when the charge velocity equ als th e spin veloc i ty. 

When we switch on interactions among spinous the above expressions ( ^.54 , 5.55, ^.58| , |5.59| ) remain formally 
unchanged if the phase shifts are suitably renormalized by Ka^ analogously to Eq. ( |5.41| ) where 5j? 6± and 
d {S -\- m). As a consequence, the exponent of the leading singularity (i.e. that corr esponding to m = 0) Xa 
exactly reproduces the exponent already obtained in the calculation of the overlap ( ( |3.42 ). Moreover, the full set of 
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expon ents Xm agrees with the term proportional to 2TrVs/£ in the finite size correction to the energy of the model 



( 5.35 ) for the particular choice of the quantum number n = 0. The rest rictio n n = is due to the fact that the 
intermediate states appearing in the rotational invariant Green function ( ^.47 ) do not have definite spin and then, 



at fixed m, the leading singularity is related to the smallest phase shift which corresponds to rt = 0. It is however 
clear how to generalize these exp ressio ns for the calculation of the spin up (or spin down) Green function: The sum 
over intermediate states in Eq. ( 5.49| ) has to be restricted to the states with the correct spin projection. Within 



bosonization it means states with the appropriate value of the qu antum number n which is in fact related to the total 
spin of the one hole intermediate state throu gh E qs. (5.£. 5.10|) . This procedure in fact reproduces the full set of 



critical exponents Xm.n which appear in Eq. ( ^.35 ) 



Finally, if we want to evaluate the singularitie s in the hole density of states we just have to integrate the leading 
momentum dependence of the spectral function ( ^.57| ): 



N{lu) = I ^A{p, u;) oc y c™ c^2X„(o)-i/2 ^ ^^^^^ (5 



where Cm are finite amplitudes, lu is measured from the bottom of the band and the exponent XmiO) coincides with 
the previously introduced critical exponent Xm evaluated at the bottom of the holon band i.e. at a total momentum 
p such that kh = p + Qm sits at the minimum of eh{kh)- In fact this region of integration in momentum space gives 
rise to the leading singularity in the density of states. Note that a divergence in the density of states occurs only if 
Xm{0) < 1/4 for some m. 



VI. SINGLE HOLE IN THE XY MODEL 



Here we analyze in some detail the dynamical properties of a single hole in the XY model on the basis of the 
Hamiltonian already introduced in Section |^ in the limit of small hole mass: J' = J «t. The aim of this study is 
to check, in a simple model, all the general features of hole motion already discussed in Section ^ and to carry out 
the quantitative evaluation of exponents and amplitudes for this system. 

The Hamiltonian of a hole of momentum p in the t — Jxy model ( [1.1| ) has been previously derived in some detail 
and reads: 



Hp — t exp 



£-1 



h.c 



(6.1) 



Here we consider a spin down hole of momentum p in a chain with even number of sites L and odd number of up 
spins (i.e. an odd number of fermions in the representation of Eq. (3.1)) N ^ 2i> + 1 corresponding to a z-axis 



magnetization fi 



+ y-. According to the discussion in Section III 



the appropriate boundary conditions of the 



fermionic problem are therefore antiperiodic, N + L being odd, and the quantization rule for the momenta is 



27r 
T 



(6.2) 



with £ — L — 1 as usual. The states which diagonaliz e th e hole kinetic term are states of given spinon momentum Q. 
In the J — !■ limit the spectrum of the Hamiltonian (|6.l|) can be obtained by diagonalizing the magnetic term in the 
subspace of fixed spinon momentum Q. This procedure gives rise to spinon Slater determinants of plane waves. The 
ground state is doubly degenerate and corresponds to a set of N occupied orbitals of momenta centered around fc = 0, 
while low energy excitations can be obtained either by changing the number of fermions (i.e. the magnetization) 
— > iV + n or by moving m fermions from the left to the right branch of the Fermi surface. The energy of these 
excited states can be easily calculated by taking the expectation value of the magnetic term on the appropriate spinon 
Slater determinant. Energy and momentum are therefore given by: 



Ep = 2tcos(p + Q) - J- 



1 



E 

j— — ij'—n-\-7n 



27r, 1, 
-0+2) 



u-\-rL-\-rn 

E 

— ~~i'—n-\-m 



COS 



27r, 1. 



(6.3) 



(6.4) 



These expressions, being based on perturbation theory in J << t, are exact to 0{J/t) and to 0(1) respectively. By 
carrying out the summations and expanding up to 0{l/£) at fixed p — N/L, we find: 
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sm irp 



J cos np (n + p) + 2t cos{p + Q) + 



27r Jsin tt/j 



Q — 2Trp{m 



(m + ^'f + ^{n + pf 



27r 



-{n + p){m 



J Sinn p- 



J cos irp 



{n + p) 



(6.5) 
(6.6) 



In the thermodynamic hmit, besides the extensive magnetic contribution, the energy of the state depends on the 
holon momentum = P + Q through the form of the holon band (fc) — 2t cos . The spinon momentum is instead 
given hy Q = 2fci?(m + i) when we recognize that the spinon Fermi momentum is just kp = t^P- The 0{l/£) size 



corrections to momentum and energy can be also compared to the general expressions (5.26,5.35) if we recall that 
the spinon velocity is Vs = J sin kp = Jsimrp and K^^ = 1 in the XY model. Equations ( |6.5D "and ( |6.6| ) correctly 
reproduce the predicted structure of the size corrections showing that the hole dynamics is described by an underlying 
conformal field theory. This comparison allows to determine, to l owest order in J/t, the values appropriate for the 
t — JxY model of the non universal parameters entering Eq. (5.35): 



for the phase shifts and 



1^(6+- S-) 



(6.7) 



27r 



(fc+ + fc^) - 



27r 



(4 



kp) = 0(J) 



(6.8) 



for the shifts of the righ t and left Fermi wavevectors. 

According to Section V C the phase shifts (|6.7| ) com pletely determine the singularities of the Green function. As 
a check, let us explicitly evaluate the overlap C ( |5.36| ) in this model. In fermionic representation this amounts to 
calculate the overlap between the ground state of Hamiltonian (6.1) and the ground state of the XY model on a 
L-site ring with the same number of fermions iV = 2i^ + 1. In ferniion representation the latter state is just a Slater 
determinant of plane waves with momentum quantization corresponding to periodic boundary conditions on a L-site 
ring: 



27r 



1] 



(6.9) 



Therefore, C is simply the overlap of two Slater determinants of plane waves with different quantization rule (6.2, |6.9 ). 
Such an overlap can be calculated as the determinant of the N x N matrix of the overlaps of the two sets of plane 
waves defined on the squeezed chain: 



Hrs =< qr\ks > = i 



detH 



COS [qr/2] 



^L{L- 1) sin [{k, - qr)/2] 



(6.10) 
(6.11) 



with r and s belonging to the interval [— i^, v\. Notice that the quantization rule of the wavevectors kj can be naturally 
interpreted in terms of a momentum dependent phase shift: 



27r, 



1^ 



(6.12) 



Near the Fermi points the phase shifts are then given by: 6± = in agreement with their determination based on the 
structure of the finite size corrections to the energy ( |6.7[ ). The known treatment of the orthogonality catastrophe 
problem thus gives th e cri tical exponent in ter ms o f the phase shift at the Fermi points: 2Xo — S"^ + 5"^ which agrees 
with the expression ( ^.39| ) derived in Section |VC| . This result can be also checked numerically by evaluating the 
determinant (5.10) for fairly large system sizes. The size scaling of In C as a function of IniV is shown in Fig. 4 for two 
magnetizations (/i — and /i = 1/4) corresponding to the densities of spin up p = 1/2 and p = 3/4. The analytical 
value of the exponent Xq is also shown in figure. The exponent for the p = 1/2 case also agrees with an independent 
calculation by Penc et al. p4|. 



The next task is the evaluation of the spin down Green function (5.44) which, via Galileo transformation, takes the 
form ( 5.461 ). In the J ^ limit, the energy levels can be written as the sum of a holon part e^ip + Q) = 2tcos{p + Q), 
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which just depends on the spinon momentum Q, and a spinon term which is eigenstate of the XY Hamiltonian on 
the squeezed chain. Therefore, at fixed total momentum p, the Green function reads 



GiP, ^) = 7 E E e-"'-(f+«)- ZiR, t) (6.13) 

Q R 

where the spinon term, imphcitly depending on the total momentum p, is 

Z{R, t) =< * I e-'"'^T^{l - no) I ^' > (6.14) 

Here He and Ti respectively represent the XY spinon Hamiltonian and the translation operator of one lattice spacing 
on the squeezed ring of £ = L — 1 sites (origin excluded) while > is the ground state of the XY model on the L-site 
ring. Notice that, for the special case of the XY model at J ^ 0, the decoupling of the Green function in holon and 
spinon factors 

is exact at all distances while, in general, we expect this decoupling is valid only at low energy and long wavelength, 
i.e. for R,T ^ 1. 

The evaluation of the matrix element in Eq. ( |6.14 ) can be performed because in the spinon representation the two 



sates are Slater determinants (without phase shifts) and the unitary operator acting on them is a one body operator. 
The projection operator (1 — rip) implies that the origin is an empty site which amounts to exclude the origin in the 
evaluation of the overlap matrix. Therefore, Z{R,t) — det A{R,t) where the overlap matrix (of linear dimension 
21^ + 1) can be easily calculated by inserting a complete set of orbitals which do not place particles in the origin. A 
useful choice is the set of eigenstates of the spinon Hamiltonian i/g- corresponding to an odd number of spinous i.e. 



with momenta (6.2). The resulting form of the overlap matrix is: 



_ COs(gj2)cOs(g,/2) ^ik,R+rJcosk,) 

The sum can be performed analytically only in the static limit r = where the numerical computation of Z{R, 0) 
and of its Fourier transform Z{Q) can be pushed to fairly large system sizes. A comparison between two different 
sizes at fixed magnetization (/i = and /i = 1/4) is shown in Fig. 5a where we can see that Z{Q) does not vanish 
outside the "spinon Fermi surface" even if it is strongly suppressed. The singularity at the Fermi momentum kp = np 
appears to be present on both sides of the spinon Fermi surface. A logarithmic plot of the s ingularity of Z{Q) when 



Q approaches kp is shown in Fig. 5b where it is compared to the analytical value of Eqs. (6.7, B.24). The ratio of 



the two amplitudes on both sides of the singularity is in rather good quantitative agreement with expression (5.6C) 



Now we examine in some detail the asymptotic form of the dynamical spinon Green function of the t — Jxy model 



(S.14). As noted before, the exact calculation of Z{R, r) in a finite system reduces to the evaluation of the determinant 
of the matrix Ars{R, t) defined in Eq. ( |6.16 ). The exact calculation can be performed only numerically and the direct 



interpretation of the data is obscured by finite size effects. However, we can address the problem of the long distance 
and long time behavior of Z{R,t) by performing the asymptotic expansion of the matrix elements A^s themselves. 
This expansion can be carried out rather easily in the low density limit where, according to expectations of conformal 
field theory, it should be characterized by phase shifts (5+ = ^_ = 1/2 (6.7) . The details are reported in appendix ^. 



According to bosonization, the spinon Green function should behave as (5.53) 



^ikpR ^ — ikpR 

^ {R-VstY/\R + v,tY/^ + {R~v,tY/^{R + v,tY/^ ^^-^^^ 

where Vg = J kp and this expression is valid for \vs t\ < R. From this analysis we expect that the function Z(R, r) at 
long wavelength behaves as 

7 (J? ^ cos{kpR) 

(i? - VsT)^I^{R + VsT)^I^ 

In Fig. 6 we plot the numerically evaluated ratio | det A(i?,r) \/Zq[Rt) as a function of R — VgT for TV = 100 and 
fi = —1/4 which belongs to the low density regime. This ratio has no oscillations (meaning that the phase factor has 
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been correctly determined) and it is approximately constant over a wide range of values of R — VsT. Clearly the region 
R ~ VsT cannot be well represented by the bosonized form which would predict a spurious divergence in the Green 
function which instead is bound to have modulus less than unity. 

The analysis of this Section shows in a simple example that all the features of hole propagation in a magnetic back- 
ground derived with field theoretical formalism are contained in such a microscopic model. This detailed calculation 
supports the assumptions introduced in our general study of the long wavelength hole dynamics. More interesting 
systems can now be investigated. 



VII. BETHE ANSATZ MODELS 

In this Section we consider two Bethe ansatz solvable models, the repulsive Hubbard model 1^ and the t — J 
model at J = 2t where the finite size corrections to the energy can be found analytically, leading to a formal 
expression for the critical exponent which appears in the one hole Green function. We analyze both models at arbitrary 
magnetization so that our results can be extended to the attractive Hubbard model via the well known canonical 
transformation. 

The calculation of the finite size corrections to the ground state energy closely follows the original derivation by 
Woynarovich for the Hubbard model at finite density Here we consider a chain of L sites and L — I electrons, 
among which have spin down. The Bethe ansatz solution is characterized by two sets of rapidities: For the 
Hubbard case we have Ng — rapidities Ac and Nc ^ L — 1 rapidities kj , while for the t ~ J model (in Sutherland 
representation) we have Ns — + 1 rapidities Va and Nc = 1 rapidity wq. These rapidities are related to the 
quantum numbers Ja and Ij respectively through the Bethe ansatz equations. In the ground state, the two sets of 
quantum numbers define two compact distributions bounded by and respectively. The explicit expressions of 
J* and are : 



J+ -J- = N, 
I+-I- = 



1+ 



J- = 2D^ 
I' = 2Dr 



(7.1) 



where Ng and Nc have been previously defined while the "centers" of the distributions Ds and Dc specify the spin 
and charge state respectively. In the ground state, at fixed total momentum p, is the smallest integer (or half 
integer) compatible with the quantization rules for J^. Instead, = Ih + L/2 for the Hubbard model and Dc = Ih 
for the t — J, where Ih defines the position of the hole in the distribution of charge rapidities and is related to the 



momentum of the holon. Following Woynarovich, we define the four "densities" as the L — > oo limit of ly^ 



and Sc(s) = Dc{s}/L, so that D, 



c(s) 



5c(s)L and Nc(s) 



is) 



Nc(s)/L 



'Ms) 



are finite for L ^ oo. If a hole of momentum p and spin 



down is created, the total momentum of the state is given by: 



2tt 

-Po + -r [{Ns 



VsL){D, - SsL) + (Nc - i^cL)[Dc - 5cL)] 



(7.2) 



where — po = 'i'T^iDc^c + DsNg)/L is the momentum in the thermodynamic limit. Analogously, the size corrections 
to the ground state energy E for the single hole problem can be expressed in terms of the above defined quantities, 
the charge (spin) velocity Wc(s), the elements of the dressed charge matrix and of the additional matrix Zij as: 



L[E - Leoo) = -77«s + '2i^VcXc + 2t:VsX 
6 



(7.3) 



Xc = -{Nc-UcL) 



{Dc - 5cL) + i^2{D, - 5,L) - AZ21 



{N,-VsL)-^^2{Nc-VcL) 
262 



X 



{Ds~5sL)~ Zi2{Nc-ycL) 



{N,-v,L)~^^2{Nc-VcL) 
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where Coo is the ground state energy per site of the model at half filling and Nc — VcL — —1 for the Hubbard case 
(where Nc = L — 1) and Nc — VcL = 1 for the t — J model (where Nc — 1). The diagonal element of the dressed 
charge matrix is simply related to the correlation exponent of the Heisenberg model introduced in Section ^ by the 
well known expression = ^22 llll • This general formula is valid both in the Hubbard and in the t — J model. The 
basic steps for the formal derivation of Eq. (7.3) together with the precise definitions of the quantities appearing in 
it are reported in Appendix ^ Here we only stress the decoupling of the charge and the spin terms in the finite size 
corrections, in agreement with our starting assumptions: The long wavelength Hamiltonian is the sum of a charge 
part which does not give any singularity in the correlation functions and a spin part which instead gives rise to critical 
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exponents. Most impo rtant ly, this rather compHcate, exact, expression (|7.3|) perfectly matches the predictions of the 



bosonization meth od ( 5.35 ) 



From Eq. (7.2) we find that the size scaling of the total momentum is characterized by the am plit ude ( 5.33| ) 



a = —{Ns — VsL){Ds — SsL) — (Nc — i'cL){Dc — 6cL) while from the magnetic contribution in Eq. (7.3) we easily 



identify the quantum numbers (m, n) a nd the phase shifts (<5, 7). The other parameters appearing in Eq. (5.35) then 



follow from the charge part of Eq. (O) leading to 



a^-{Ns~ v,L){D, ~ S,L) - (TV, - i/ci)(i?c - ScL) 

m = Ds ~ - 
2 

n = {Ns - VgL) 

w+ = 

/^Z2i{N,-v,L) 
w_ 

?22 



2Zi2{Nc-VcL) (7.4) 



where we used that in the ground state Ds = 1/2 and then 5s = 0. Note the vanishing of lu+ which occurs in all 
the models we have examined and is probably related to a sort of Luttinger theorem which forces the volume of the 
spinon Fermi surface to be unaffected by interactions. This exact correspondence between bosonization and Bethe 
ansatz demonstrates the validity of our approach in the Hubbard and t — J models and allows for the analytical 
determination of the phase shifts 5± governing the sin gularities of the one hole Green function. 



In general, the coefficients appearing in Eq. (7.3) are non universal quantities which depend on the coupling 
constants of the model as well as the average magnetization per site /i and the holon momentum. Therefore, we 
expect that the exponent X is a function of all the parameters which define the hole Hamiltonian, including the total 
momentum p. However, some special but important exception must be mentioned. At zero magnetic field the spin 
Hamiltonian possesses the additional SU (2) symmetry both in the Hubbard and t — J model. At this particular point, 
the exponent X is universal. In fact, the elements of the matrices ^ and Z acquire analytic values independent of the 
total momentum of the state and of the coupling constants: ^22 = l/v^, ^12 = 1/2 and Z12 = 0. Another simple case 
is the U ^ 00 limit of the Hubbard model where Z12 = and the dressed charge matrix can be expressed in terms of 
the magnetization per site — {N^ — Ni)/2L: ^12 = 1/2 — while Kg- = £,22 as a function of fj, is shown in Fig. 1 of 
Ref. [|2|. Also in this limit the exponents do not depend on the total momentum of the state. Finally, when the hole 
sits at the bottom of the band, i.e. if we are at the one hole ground state, the equations simplify because the holon 
momentum is always kh ~ it and so Vc ~ 0. Again Z12 = and our expression for the finite size corrections coincides 
with the zero doping limit of the known form valid at finite density. This proves the continuity between the physics 
of the single hole problem and that of finite doping in ID. 

At arbitrary magnetization and momentum no analytical expression for the phase shifts is available. However, the 
integral equations reported in the A ppen dix p| can be solved numerically. Few examples are reported in Fig. 7 where 
we show the critical exponents Xm ( 5.42| ) as a function of the momentum for several magnetizations in the Hubbard 



model. The loci uj = ^h{p + Qm) in the (p, uj) plane where the occurrence of the divergence is predicted by our theory 
are instead shown in Fig. 8 for a couple of choices of the parameters. The presence of 2kF zero energy excitations in 
the spinon spectrum gives rise to a remarkable symmetry property of these curves: At a given energy w, if a singularity 
occurs at momentum p it will also show up, with possibly different exponent, at momentum —p + 2kp. 

As a check on the theory of Section ^ we have numerically evaluated the overlap of Eq. ( p. 36 ) by Lanczos 



diagonalization in the t — J model at J = 2t in chains up to 32 sites at magnetization /i = ±1/4. The results are 
well fitted by a power law behavior in L ( |5.38D with an exponent X which clearly depends on the momentum of the 
hole, in agreement wit h th e bosonization analysis. A comparison between the numerically determined exponent and 
the prediction of Eq. (|7.3| ) is shown in Fig. 9 (a). An analogous computation has also been performed in the < — J 
model at the generic non integrable point J — t and /i = 1/4. The numerical results are also shown in figure 10 (b) 
but in this case the comparison with the analytical predictions based on the Bethe Ansatz solution is not available. 
However, a quite similar momentum dependence of the overlap exponent X emerges, showing that the above features 
are not special to the exactly integrable points. 
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VIII. THE SHAPE OF THE SPECTRAL FUNCTION 



In this Section wc finally discuss the global shape of the spectral function of one hole in a correlated background. 
In fact, the bosonization method developed in Section ^ only concerns the long wavelength, critical properties of the 
Green function and gives no information on its short wavelength features. The exact calculation of a Green function 
in an interacting system, at all lengthscales, has been achieved only recently |^ by use of extremely sophisticated 
methods for the Calogero-Sutherland model |^ while no results are available in other interacting systems. Our task 
is twofold: to understand the physical nature of the low lying excitations which contribute to the spectral function 
and to develop a useful numerical method for the approximate determination of the spectral function A{p, lu) which 
can be applied to generic one dimensional models. The purpose is to overcome the severe finite size effects present in 
Lanczos diagonalizations js^js^ without resorting to the delicate extrapolations of simulation data |Q necessary for 
the computation of dynamical correlation functions. 

Here, we formulate a simple approximation for the one hole Green function which captures most of the features 
of the exact result and can be usefully applied to interesting correlated models, in one dimension, like Hubbard. A 
similar approximation has been discussed by Penc et al. |23|| in the strong coupling limit of the model. For clarity, 
we introduce the method in the framework of the previously discussed t — Jxy Hamiltonian in the J — *■ limit (see 



Section VI ) leaving the study of the Hubbard model as a final example. 



The starting, exact, expression f or th e ID Green function of the t — Jxy model is Eq. (5.13) where the spinon 



function Z{R,t) is defined by Eq. (6.14). The spinon term Z in fact contains all the interesting correlation effects as 
previously pointed out. However, its direct evaluation proved rather hard even in the simple case of the XY model 
where the exact ground state is a Slater determinant, while the more realistic c ase of a Hubbard model cannot be 



tackled by these methods. As an approximate way to evaluate the matrix element ( 6.14 ) we can assume that the most 
relevant contribution to the intermediate states comes from the single spinon states, i.e. from the exact eigenstates 
of H„ in the squeezed chain with only one spinon excitation. In fact, being an antifcrromagnetic Hamiltonian 
defined in an odd chain with periodic boundary conditions, it gives rise to a frustrated problem. Then, its ground state 
contains a free spinon and it is rather natural to assume that a set of low energy states can be built by giving a finite 
momentum Q to such a "quasiparticle" . These single spinon states \Q > can be therefore labeled by the momentum 
Q, which lies outside the spinon Fermi surface, and have energy ea-{Q) given by the spinon band: ea-{Q) = JcosQ for 
the XY model, where \Q \ > kp and kp = np is the spinon Fermi momentum. The single spinon approximation to 
the function Z{R,t) therefore reads: 

ZiR,r)^j e-'^"(^)"e'Q-«| < Q|(l-no)|* > P = i ^ e-''^^^^^ e'^'^ Z{Q) (8.1) 

\Q\>kF |Q|>fcF 

where | ^' > is the ground state of the undoped model and we have introduced the spinon function Z{Q) 

Z{Q) = e^"'^'^ < * I Tf {l - no) I * > (8.2) 

R 

This approximation, therefore entirely resides in the assumption that only the single spinon intermediate states give 
a finite contribution to Z(R,t). As a consequence, the time dependence of Z{R,t) is greatly simplified but still not 
trivial due to the complex structure which can be present in Z{Q). 

The value of this approximation is that the time dependence of Z{R, r) is given analytically in terms of the known 
spinon excitation spectrum of the model and that only a small 0{t} number of matrix elements is necessary for the 
evaluation of the full Z{Q). The exact calculation in fact would require the i nser tion of a complete set of intermediate 
states leading to an exponentially large number of terms in the sum of Eq. ( |8.lD . A consequence of the single spinon 
approximation is the presence of a sharp Fermi surface in the function Z{Q) which is in fact predicted to vanish 
identically for \Q\ < kp. This property is well satisfied in the quasiparticle weight for the J ^ limit of the t — J 
model JTl| , p3t even if both the numerical calculation of Z(Q) in the t— Jxy model (see Fig. 5a) and its density matrix 
renormalization group evaluation in the strong coupling limit of the Hubbard model show that this is not an exact 
feature of one dimensional systems [Q . 

We now proceed to the evaluation of Z{Q) in single spinon approximation for the exactly soluble t — Jxy model, in 
order to provide a check on the quality of this approximation. The completeness condition of the intermediate states 
which is obeyed by the exact quasiparticle weight reads: 

Yz{Q)^i{l-p) (8.3) 
Q 
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In order to compute Z{Q) in the t — Jxy model on a L-site ring and 2i/ + 1 spinons, we first recall t hat t he ground 



state of the XY model is a Slater determinant of 2i/ + 1 plane waves with momentum quantization (6^) and then 
the occupied orbitals are | qs > with s_G [— i', i']. The intermediate states are also Slater determinants of plane waves 
but with different quantization rule (6^). If the spinon momentum is Qr — 5)' ^^'^^ single particles 



intermediate states | kj > fill the spinon Fermi sea j € [— z^, u — \\ while the remaining spinon is placed outside this 
interval, at j = r. As noted before, the overlap between two Slater determinants is just the determinant of the matrix 
of the overlaps Bsj —< qs\kj > and the quasiparticle weight Z[Qr) is the modulus square of such a determinant. 
The numerical evaluation of the completeness sum rule ( p.3[ ) in single spinon approximation is plotted in Fig. 10 for 
two magnetizations. The data show that the sum rule is violated in the thermodynamic limit and therefore the single 
spinon states do not represent a complete set of intermediate states as expected. However, Fig. 10 also shows that 
the breakdown of the sum rule is very small and shows up at considerably large system size. This approximation 
accurately reproduces the short wavelength properties of the model while fails in catching the long wavelength features 
(i.e. critical exponents and amplitudes) which we already discussed by use of bosonization methods in Section 
Therefore, we expect that single spinon approximation can be successfully applied to the study of the global shape 
of the one hole spectral function in one dimensional models. As an example we now briefly discuss the case of the 
Hubbard model. 

Analo gous ly to the t — Jxy model, the Green functi on is expressed in terms of the dynamical quasiparticle weight 
by Eq. ( |6.13 ) which is approximately evaluated as in (8.1): 



G(p, r) = i ^ e-^^l.{p+Q)r^-^cAQ)r\ < Q | ^^^^ | ^ > |2 (g^4) 

Q 

The remaining problem is to compute the function Z{Q) — (. \ < Q | Cp^j^ | Vl/ > p for the Hubbard model. The 
numerical calculation can be carried out by Lanczos diagonalizations by exploiting the negligible size dependence of 
this quantity already verified in the previous examples. The ground state > of the half filled Hubbard model at 
magnetization = — i + ^x"'" ''^^ numerically obtained in chains up to L = 16 sites. The remaining problem is to 
select the one spinon states \Q > which contribute to Z{Q). The procedure we have adopted takes advantage of the 
continuity of the one spinon states between the weak and the strong coupling limit. At [/ — s- cx) the one spinon state 
\Q > is in one to one correspondence with the ground state of momentum Q of the Heisenberg chain in a ^ = L — 1 
site ring. This follows from the factorization property of the Hubbard eigenfunctions discussed in Refs. |p^ , |ll| . These 
spinon states can be therefore identified by performing Lanczos diagonalization on the £ site Heisenberg model in the 
symmetry subspace of total momentum Q. Now, having determined the one spinon states at U = 00, we adiabatically 
lower the interaction parameter U (in practice this is quite an easy procedure within Lanczos method) following the 
"evolution" of the eigenstate of spinon momentum Q from U = 00 down to the desired value of U. In this way, 
starting from an exact eigenstate of the Heisenberg model we first find an eigenstate of the Hubbard model at J7 = 00 
and then a sequence of eigenstates of the £ site Hubbard model corresponding to smaller and smaller interaction 
parameters. This procedure has been devised because the spinon momentum Q is a good quantum number only at 
U = 00 and then we need a method to select the exact one spinon states out of the full set of eigenstates of the 
Hubbard Hamiltonian. As usual, a simple check on the validity of t he s ingle spinon approximation comes from the 
completeness condition of the intermediate states, analogous to Eq. ( |8.3| ), which now reads 

i Z{Q) -< * I cl^cp,i I * >EE ip) (8.5) 
Q 

i.e. the momentum distribution of the spin down electrons at half filling. The amount of violation of this sum rule 
quantifies the weight of all the other states in the Hilbert space which have been neglected in our approximation. 

In Figs. 11-12 we show results from Lanczos diagonalization in chains up to 16 sites. The function Z{Q) for two 
values of the momentum p of the hole and of the magnetization /i is reported. In all the cases we have considered the 



completeness condition (8.5) is very well satisfied showing that one spinon states account for more than the 98% of 
the full Hilbert space of intermediate states js^] . This allows to reconstruct the full spectral function for the Hubbard 
model. On the other hand, the available data also show weak size dependence suggesting that finite size effects are 
not relevant at high energies. A plot of the predicted spectral function is shown as a function of uj in Fig. 13 for a 
typical, intermediate coupling {U = At) and a two hole momenta p. 

IX. SUMMARY AND DISCUSSION 

In this work we have analyzed in some detail the dynamical properties of a hole in an antiferromagnet. Due to 
the mapping between the attractive Hubbard model at arbitrary density and the half filled repulsive Hubbard model 
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in a magnetic field, our analysis directly applies also to the more general case of hole propagation in correlated one 
dimensional models with a gap either in the charge or in the spin spectrum. Most of the known quasi one dimensional 
materials in fact belong to these classes and then the present study may be helpful in the interpretation of the available 
photoemission spectra of quasi Id systems First we found the exact spectral function of a single hole in the Ising 
model which is characterized by a gap both in the charge and in the spin channel. As a result, the quasiparticle weight 
is finite also in d = 1 and A{p,uj) has a 6 contribution. However, due to the absence of spin fluctuations the hole 
dispersion relation is flat and the hole cannot propagate. Then, we focused on the singularities of the hole spectral 
function which occur because of the presence of other gapless degrees of freedom (spinous) . This spinon gas behaves 
as a Luttinger liquid which gives rise to the typical critical exponents of one dimensional physics that show up in 
the hole dynamical properties. Among the results we have obtained, we like to stress few general features of A{p, lo) 
which characterize hole propagation in one dimension: 

• 1) The main singularities in the spectral function occur along lines in the (p, w) plane with dispersion relation 
determined by the form of the holon band, uj{p) = th{p^ kp), while the spinon excitation induced by the hole is 
created at the Fermi points ikp- This is a consequence of spin charge decoupling which occurs in one dimension 
and gives rise to divergences in A{p,Lu) also above the bottom of the band. In this case, the divergence may 
occur on both sides of uj{p) with different amplitudes. The existence of two branches of singularities in the 
spectral function (see Fig. 8) can be interpreted as due to the presence of a shadow band psf . 

• 2) The singularities are characterized by critical exponents which can be explicitly calculated in integrable 
models. In the isotropic antiferromagnet the SU{2) symmetry forces the exponent to be exactly X = 1/4 for 
all microscopic Hamiltonians. Instead, when spin isotropy is broken, or when the system has a spin gap, the 
critical exponent X in general depends on the parameters of the model and on the momentum of the hole. 

• 3) Away from the SU{2) isotropic point, the tunneling density of states has either a divergence or a zero at the 
bottom of the spectrum according to the value of the critical exponent X, i.e. according to the parameters of 
the model. Remarkably, the density of states shows at most weak logarithmic singularities in the isotropic case. 

• 4) At the bottom of the band the critical exponent coincides with the known exponent characterizing the spectral 
weight in doped systems when the zero doping limit is taken. This provides a demonstration of the continuity 
of the physical behavior of the degrees of freedom which do not develop a gap in the excitation spectrum as 
doping vanishes. 

All these features should be experimentally detectable in quasi one dimensional materials. 

A problem posed by this analysis concerns the relationship between our results and the zero doping limit of the 
generally accepted Luttinger liquid picture of the Hubbard or t — J model. In fact, the continuity between the 
single hole and the low doping physics apparently breaks down when the hole momentum does not coincide with 
the Fermi momentum of the doped model, i.e. when we are above the bottom of the holon band. In this case, 
standard bosonization methods would predict singularities in the spectral function with momentum independent 
critical exponent uniquely determined by the physics at the Fermi points p7[ | while the accurate analysis of the single 
hole problem reveals the presence of momentum dependent critical exponents when we move away from the Fermi 
level. 

From a methodological point of view, this study of the single hole motion demonstrates a close relationship between 
the physics of hole motion and the single impurity problem in Luttinger liquids. This mapping is provided, at 
strong coupling, by the Galileo transformation which allows to eliminate the hole degree of freedom in favor of a 
non translationally invariant spin system. The recoil of the hole, embodied in the hole kinetic contribution of the 
effective spin Hamiltonian iJp, cuts off the backward scattering terms in the impurity problem and generates effective 
boundary conditions which allow the propagation of spinous through the impurity site. This idea can be extended 
to higher dimensions. The Galileo transformation, in fact, is not restricted to d = 1 and the single hole problem can 
be always mapped to a pure spin Hamiltonian. Generalizing what we found here, it is tempting to assume that also 
in d > 1 the hole acts as an effective boundary condition placed at the origin of the d-dimensional spin lattice. The 
emerging picture resembles and generalizes that of the dipolar distortion proposed by Shraiman and Siggia |Q based 
on the semiclassical treatment of a particular choice of boundary condition. More work along these directions may 
eventually clarify the properties of hole motion in a correlated background in arbitrary spatial dimensions. 
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APPENDIX A: FORMAL CALCULATION OF THE OVERLAP IN A LUTTINGER LIQUID 



In this appendix we give a formal derivation of the overlap 

0±((5±)=<d±|c±> (A.l) 

between two free particle states 

< d±\= II di\0 >, <c±\=l[cl\0> (A.2) 

on a given branch of a Luttinger liquid, here identified by (±). In order to simplify the notation we drop the label ± 
from the operators c and d whenever it does not lead to ambiguities. The operators d are defined with skew boundary 
conditions (i.e. with non vanishing phase shifts) while the operators c correspond to periodic boundary conditions: 

£ 
2 

_ i_ 

2 

£ 
2 

di = ^ J e*T(«+A')-^^(a;) (A.3) 

2 

where V'±(a;) identifies the fermion field in the right (+) or left (— ) branch. The relationship between the operators 
d and c is easily found using canonical anticommutation rules for the fields tIj±: 



oc 

4 



n— — oo 



J2 ss{n--m)dl (A.4) 
where 

ssin) = (A.5) 
■7T(n + 0) 

The overlap 0±{S) is given by the determinant of the overlap matrix Dn^m{5) = ss{n — m) with the restriction on the 
allowed indices n, m > for the left branch and n,m <Q for the right one which selects the occupied orbitals. From 
these definitions we get the symmetry property: 

0+{5) = 0^{-5) (A.6) 

which follows from the transformation rule of the matrix Dn^m{S) under the mapping (n, m) (— n, — m) which 
changes the left into the right branch. A further property of the determinants 0±{S) derives from the definition of 
the matrix D{6): -Dn,m(<5) ~ ss(n — m) ~ S-s{m — n) = D„i.n{—S) which gives: 

0±{S)^0±{-S) (A.7) 

as the determinant of a matrix is equal to that of its transpose. In the continuum limit the determinant 0± is not well 
defined and a cut-off procedure is required before evaluating the overlap. Some care should be taken in the explicit 
definition of the cut-off. In fact, by restricting the matrix indices (n,m) to a finite interval we would effectively 
introduce an unphysical doubling of the Fermi surface. Instead, let us consider a system with both left and right 
branches with the same finite phase shifts (5+ = (5_ = 6, and finite but large number of particles N symmetrically 
distributed in the positive and negative branch. In this limit it is clear that the left and right branch decouple and 
the total overlap 0{5) is given by the product of the two left and right component 0{6) ^ 0+{S)0-{S). The overlap 
0{S) can be exactly computed in a finite lattice with a given (large) number of particles N and is formally given 
by the previously introduced determinant with matrix indices belonging to the interval [1, A^]: 

Oi5)=aet-^^^=A,N--^' (A.8) 
7r(n — m + 0) 



where Ag is a finite numerical constant. Then, by use of the relations (A. 6, A. 7) and the previous result we finally get: 

I O+iS) P = I 0-i6) \' = 0(6) = AsN-'' (A.9) 
which gives | 0±{S) \ oc £^^^^ where we expressed the number of particles as a fraction p of the number of sites £. 
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APPENDIX B: THE DYNAMICAL QUASIPARTICLE WEIGHT 



In Eq. (5.52) we showed how the dynamical quasiparticle weight Z{x,t) sphts into the product of contributions 
coming from the right and left fermion branches. Now we have to evaluate the generic matrix element appearing in 
the formal expression of Z(x,t) for a non interacting Fermi gas characterized by given phase shifts 5± at the two 
Fermi points: 



Zi{x) c± le"^^*-'^*) |c± > 



(B.l) 



where the states \c± > are defined in appendix A an d P± represents the total momentum operator for the fermions 
in the right + or left — branch, defined in Eq. ( |3.51 ), with fermionic operators '0+ ^^nd ip- obeying skew boundary 



conditions (5.4). Here (5-|_ is the minimum momentum of the right branch of d electrons (i.e. electrons with a defined 
value of the momentum P ) , as excitations in the right d branch can only increase the momentum by ^ with positive 
integer j. Conversely, Q_ is the maximum allowed momentum in the left branch of d electrons, as excitations in the 
left d branch can only decrease the momentum by 2ii with negative integer j. Though in a Luttinger liquid Q± 
are infinite constants, the functions Zj^(x) and Z-(x) are finite also in the continuum limit and have an important 
property, referred in the following as the Fourier positiv ity ( negativity). By inserting a complete set of the mentioned 
excitations with definite momentum in the RHS of Eq.(B.l) we get 



n>0 j 



where j in | dj^n > labels all the possible excited states | dj_„ > with momentum (5+ 



(B.2) 



Such a spectral decompo- 



sition of Z^j^ implies that the Fourier coefficients (n) = 



' Z^{x)dx are non vanishing (and positive definite) 



only for n > (Fourier positivity). Analogously Z^{n) > only for n < (Fourier negativity ). 
In terms of the d operators (|A.3|) the momentum on each branch P± is diagonal and reads: 



(B.3) 



where k, 
transformation n 



^{n + S) and kp — ±kp. Using the same par ticle hole transformation discussed in appendix A, (i.e. the 



-n) we obtain the analog of Eq.( A.6) 



Ziix) - Z-'{-x) 



(B.4) 



Having discussed the general symmetry properties of Zj.{x), we now turn to the explicit evaluation of the function 



Z\x) =< ^-le'-f^^l* > 



(B.5) 



in a free Fermi gas of TV particles, where | > is the ground state with periodic boundary conditions while the 
momentum operator P refers to a system with skew boundary conditions. Following appendix we carry out the 
calculation for a model of fermions with constant phase shift 6 throughout the Brillouin zone and then we relate 
the result to Z^{x) by use of the symmetry properties previously discussed. In particular, if we keep only the most 
relevant singularity, the factorization property proved in Eq. (5.52) gives: 



Z\x) = e'^"" Zlix) Zi{x) 



Z%{x)Zl\-x) 



(B.l 



where use has been made of the symmetry (B.4) and Q = (5+ + Q_ is the reference momentum of the intermediate 
states. 

In a free Fermi gas, both the ground state | 5* > and the intermediate states are Slater determinants built with 



the different fermionic operators c„ and d„ respectively (see Eqs. A. 2 and A.3). The translation operator exp(iPa;) 
with P — P+ + P_ ( |B.3| ) is a one body unitary operator which maps the Slater determinant with plane waves single 
particle orbitals 



0„(O = ^e^ 



(B.7) 



into another Slater determinant with orbitals 
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■j (x+r) 



where ssin) is defined in (A. 5). The overlap between these two Slater determinants is just the determinant of the 
N X N matrix M„,„ of the overlaps between the occupied orbitals: 



(B.9) 



giving explicitly 



-1) 



n — m giTTf^ 



sin ttS z" — 



1 + me 



TT 72 — m 

c siuTT^ 



for n ^ m 



for n = m 



(B.IO) 



where we have introduced the phase factor z = with tt = The matrix M„m can be written as the product 
three diagonal matrices (with diagonal elements (—1)"", and (—1)™) which contribute to the determinant with a 
phase factor and the Tocplitz matrix T„_m given by: 



de 

2^ 



C{e) = 1 + Q{u - 0){e'^''^ - 1) 



(B.ll) 



where Q is the step function. The leading singularity of the determinant of T might be extracted by means of Szego's 
theorem [p5| which would give: 



det T = exp Ngo + ^ ng_„g„ 

n>l 



(B.12) 



with 



<7n = £^ge-™« \nC{e) 



By performing the Fourier transform, gn is simply evaluated: 



go = iS u 

g„ = f(l-e-™") forn^O 



(B.13) 



(B.14) 



However, this theorem holds only for continuous C{9) and does not apply directly to our case. However one can follow 
the same regularization applied in an analogous calculation for the Ising model pst by noting that the Hilbert matrix 
of elements i?„_,„ defined by: 



Hr, 



sin nS 

n(n + 6) 



(B.15) 



has the same kind of singularity shown by T. In fact, the coefficients g„ which characterize H are just given by 
gn = for n 7^ and g^ — 0. Therefore, following Ref. we can apply Szego's theorem to the ratio of 

determinants 



detT 



^ exp <^ Ngo + ^ : 



n>l 



2,52 



9-n9n 



(B.16) 



Finally, the determinant of a Hilbert matrix can be analytically evaluated, as shown in Eq. (A.S), giving, to leading 
order, 



Z^{x) = det T = e'Q'' g+{x) g+{-x) 



(B.17) 
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where the overall phase factor depends on the reference momentum of the intermediate states Q = J^n (f*^'" ^" ^ ^f) 
and 



exp 



n>l 



= As[Nil-e'-)] 

As is the same numerical constant appearing in Eq. (|A.^) and is given explicitly by 



In^A 



C 0.5772 . . . being the Euler constant. 



-isni+c) 



(B.18) 
(B.19) 

(B.20) 



The Fourier positiv ity prope rty of g~^{x) can be easily proved by expanding the exponential in Eq. ( B.18 ). By 
comparing Eq. (B.17) and Eq.(B.6) we find: 



g+{x)/Zi{x)=g+{^x)/Z-'{^x) 



(B.21) 



where the left hand side has the Fourier positivity property (as the ratio of two function satisfying this property also 
satisfies the Fourier positivity property) , while the right hand side has the Fourier negativity property. Therefore both 
terms of this equation have to be constant in x implying that Z^{x) oc o+(x). The overall proportionality constant 
can be determined by noting that the n — Fourier coe fficient in Eq. ( [B.2D coincides with the square of the overlap 
|0+(5)P which has been explicitly calculated in Eq. (|A.9| ). On the other hand, the n — Fourier component of g^{x) 
can be read off from Eq. ( B.18[ ) yielding Z^{x) — g^{x). 
In the thermodynamic limit 



oo at fixed x and fixed density of fcrmions p = jj) Eq. (B.IE) simplifies: 



Zlix) ^ As{2Trp)-^ i-ix + e)- 



(B.22) 



with e ~ 0{x'^ / 1^) is a vanishingiy small positive term which defines the branch cut for the non integral exponentiation 
of the complex number {—ix + e). Taking the Fourier transform, we get (for < (5 < 1): 



(B.23) 



The asymptotic evaluation of the two functions Z^{x) for a Fermi gas with arbitrary phase shift also provides 
the leading singularity for Z^''^'^"^\x) = zf'''^"^\x) entering the spectral function of the interacting model with 



different phase shifts on the two branches 6±. In particular, for t = 0, the Fourier transform of Z{x, 0) in Eq. (5.52) 
splits into the sum of terms with different phase shifts 5+ = (5+ + m and S- = 6- + m each given by 



Z{Q + k) = B sininSDk 
Z{Q + k) = Bsm{TrP){-kY 



for k> 
for k < 



(B.24) 



if Q is the average spinon momentum. The prefactor can be also calculated: B = 2(27rp)-(^++^-) Ag^Ag_T{l-dl-S^_). 
In this case, the Fourier transform is nonvanishing both for positive and negative k and shows singularities with the 
same exponent for k ^ 0^. The amplitude ratio tends to a number which, in the free Fermi gas, just depends on the 
phase shifts: 



— lim , „ 
A^ k-.o+ Z{Q - k) 



Z{Q + k) _ sin(7r(52 ) 



(B.25) 



sin(7r5^) 

The asymmetry between the two sides of the singularity is therefore enhanced when one of the phase shifts gets small. 



APPENDIX C: THE LONG WAVELENGTH EXPANSION 



The first task is to evaluate the summation present in Eq. (3.16) at large R. For simplicity we restrict the analysis 
to the low density regime. It is easy to verify that the most relevant contributions to the sum are those around kj 0. 
For r ^ s this gives the approximate form: 
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where the function f(q) is given by: 



Ars - J^^—-^ [f{qr) - f{qs)] (C.l) 



r rlh pi{kR-Tjk^/2) 

Standard asymptotic expansion leads to 

/(g) = e^(9fl-r V/2) , ^ (C 3) 

J-00 27r X 

where a = [z/{l + gz)^]/(2i?) and z = —Jt/R. In the i? — > cx) limit at z =const, a — > and the problem simplifies. 
Notice that this last limiting procedure requires that a is regular and therefore that 1 + is always positive. This 
implies \z\ < l/kp (where kp is the Fermi momentum of the spinon) which is always satisfied at low density. In 
general this inequality leads to R > VsT where Vs = Jkp is the spinon velocity at low density. In the following we will 
consider only this regime. In this case, a can be set equal to zero in the integral leading to the final expression: 

f[q)=' e'i'iR-JiV^) (C.4) 



which inserted into (C.l) gives 



7r(r — s) 



^i{xr+yr ) _ pi{xs+ys ) 



where r and s run over the occupied spinless fermion orbitals [— i^, u] and x — 2ttR/ L, y — —2tt^Jt/LP'. The diagonal 
elements require a separate analysis which gives 



4 ^ piJr i{xr+yr'^) 



^ a; + 2yr 



(C.6) 



Now it is convenient to express the matrix A^g as the product of a real matrix B^s and other diagonal matrices 
Drr = e*(^'~+s"'')/2. In fact, A = DBD where 

2 sin [Ijxjr - s) + yjr^ - s^))] 
Brs ^ -, ^ r / s 

TT[r — S) 

B„^l^ (C.7) 

TT 

The asymptotic form of the A^-spinon Green function Z{R, r) is then given, besides a global phase factor coming from 
the determinant of 13, by the determinant of the real N x N matrix Brs ■ For convenience, we fix a finite ratio N/ L — p 
(even if the expressions previously derived are exact only for p — *■ 0) and we perform the numerical computation of 
dot B for several values of R and r. The results are reported in Fig. 6 

APPENDIX D: DERIVATION OF THE FINITE SIZE CORRECTIONS 



In this Appendix we sketch the derivation of the general formula (7.3) for the finite size corrections in the Bethe 
nsatz soluble models: Hubbard a 
ork on the finite size corrections 
models will be treated separately. 



ansatz soluble models: Hubbard and t — J at J = 2t. We closely follow theprocedure detailed by Woynarovich in his 
work on the finite size corrections of the Hubbard model at finite doping [g^ . For sake of clarity, Hubbard and t — J 



1. Hubbard Model 



We consider the Hubbard Hamiltonian on a L site chain at fixed chemical potential fi and magnetic field h (in units 
of the hopping amplitude t) : 
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L L L ^ L 

i—1 (7 i—1 i—1 i—1 

where a =t, [ is the electron spin index. The Bethe ansatz equations for the Hubbard chain read: 



sin kj — \f3 

fi=i 



Lkj — 27r/j + 2 arctan ( 4 



U 

2 arctan U ^" \ = 2nJa + X] ^ arctan (2 ^ '^^ j (D.2) 

j=l ^ ^ /3=1 ^ ^ 

where A'c = L — 1 and — iV^ . The quantum numbers Ij and Jq, are integers or half odd integers depending on 
the parities of Nc and Ng-. Ij = iVs/2mod(l), Jq, = {Nc + Ng + l)/2mod(l). The existence of a solution to these 
equations requires that each set Ij and Jq consists of mutually different quantum numbers. Therefore, the distribution 
Ij is uniquely defined by the position of the hole I^- Due to the periodicity of the Bethe ansatz equations by the 
substitution Ij Ij + L and kj — > kj + 2ti we can always assume that the i — 1 quantum numbers Ij fill the range 
[J/j + 1, //i + L — 1]. The important low energy real solutions are found if the Ns quantum numbers J a are chosen 
as contiguous integers (or half odd integers). We denote as Jmin and Jmax the minimum and maximum value of the 
distribution respectively. The ground state corresponds to the most symmetrical distribution around zero compatible 
with the quantization rules. It is useful to introduce the additional quantities I^ = I h — 1/2, = Ih + L — 1/2, 



J = Jmin — 1/2 and J+ — Jmax + 1/2 which, by definition, satisfy the relations ([7.l[). Finally, the total energy is 
expressed in terms of the rapidities kj and Aq. by 

E = -2Y^ cos kj + fiNc + h{Ns - Nc/2) (D.3) 
i 



The Bethe ansatz equations (D.2) can be written as Zc{kj) = Ij/ L and Zs{\a) = Ja/ L in terms of the functions 
Zc{k) and Zs{\) defined by Eqs. (2.6) of Ref. [^. Following this work we also introduce the boundaries k^ and A^ of 
the rapidity distributions defined by Zc{k^) = I^ / L and Zs(A^) = J^/L. The distribution function for the rapidity 



Pc{k) is defined as the derivative of Zc{k) with respect to k and, due to Eq. (7.1), satisfies 

Ik 



^N'%c(fc)-^' Pc(fc)j-^^' 2 arctan(4A/C/)p,(A) = ^ - 1 



(D.4) 



while the analogous distribution fimction Ps{X) obeys Eqs. (2.14) of Ref. Following the derivation of Ref. |27| ] 

we find the 0(1/L^) correction to the rapidity distributions which can be expressed in compact notation as: 

p(fc,A) =poo(fc,A)+ (D.5) 
1 f pi{k,\\k+\+) pi{k,X\k-X-) p2{k,X\k+X+) p2ik,X\k-X')' 



24L2 V Pcik+) pc(fc-) PsiX+) Ps{X-) 

where p, pi, p2 are vector functions with two components, one referring to charge and the other to spin, and satisfy 
equations (2. 16), (2. 20) of Ref. ||2^]. Notice that in our one hole case the functions pi are symmetric under the 
simultaneous interchange of fc+ — s- k~ and A+ A^ while the functions p2 are antisymmetric. This property leads 
to the cancellation of the charge contribution to the finite size corrections of the total energy which now reads 

E = Le^Xk+,k-,X+,X-) - -l-e2(fc+, A;-, A+, A-) (D.6) 

where eoo and €2 are defined by Eqs. (2.23)-(2.25) of Ref. |2^. In the thermodynamic limit k~ kg (i.e. the holon 
momentum) and k^ ^ kg + 2tt while, as usual, the spin rapidities A are centered symmetrically around the origin: 
A^ — > Aq and A^ — > —A''. By expanding eoo(fc^, , X'^ , A^) around its limiting value for infinite size we get: 
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eoo{k~^,k ,A+,A ) = £00(^0 + 27r, fco, Aq, -Ao) + 

Trvc{ko)pocc{ka) [(fc+ - fco - 27r)^ - (fc" - fco)^] + 
7TVs{Xo)poosM [(A+ - Ao)2 + (A- + Ao)2] 



(D.7) 



in terms of the charge (wc(fco)) and spin (ws(Ao)) velocities. The next step is to express the difference between A^ [fc^] 
and its asymptotic value Aq [fco] in terms of the known parameters Afc(s) a-nd i?c(s)- To this end we start from the Eqs. 
(D.4) and evaluate their derivatives with respect to fc* and A^ in the thermodynamic limit. The final expressions 
for i^c and i/g coincide with Eqs. (2.34) of Ref. in which the fc integration is extended to the full interval [0, 2n]. 
Instead, the equations for 6^ and 6s now read: 



dk+ 



dk- 

Pooc(fco 



Pooc(fco) 



27r 



(Tic(fc) 



ko 



CTic(fc) 



Pooc(fco) 



arctan(4A/C/)CTis(A) 



-An 



dk+ 



dX- 



Pooc{ko)Zii 

dSs Poociko) 



dk- 

Poos(Ao) 



2 

Poos(Ao) 
TT 

Poos(Ao) 

2 

Poos(Ao) 
TT 

dSs 



2 

27T 



o-is(A) - 



Ao 



-An 



O'ls(A) = Pooc{ko)Zi 



0'2c(fc) 



d\+ dX- 
where the additional functions 



arctan(4Ao/?7) 

fco 

0'2c(fc) - 

:7r 

arctan(4Ao/?7) 
Poo£(Ao) 



(T2c{k) 
Poos(Ao 



arctan(4A/[/)CT2s(A) = /Ooos(Ao)^2i 



o-2c(fc) 



-fco 

Poos(Ao) 



arctan(4A/[/)CT2s(A) ^ poas{Xo)Z2i 



Ao 



0-2s(A) 



-Ao 



0-2s(A) 



,(Ao)^2 



O'lc(s) 
0'2c(s) 



5p, 



9fc+ 
dpcjs) 
dX+ 



(D.9) 



satisfy the equations (2.16), (2.36) of Ref . 



From their definitions, it is easy to see that the equations satisfied by 



the elements of the dressed matrix simplify and give = 1 and ^21 = while 

r-An 



/An J\/ 
-lif2(A-A')62(A') 
-An 



62 (fc) 



-A 
dX 



-A( 



27r 



ifi(sinfc-A)e22(A) 



(D.IO) 



where the kernels Ki{x) and ^^2(2;) are defined by: Ki{x) = 8U/{U'^ + 16^^) and K2{x) = AU/{U'^ + Ax'^). Some of 
the elements of the additional matrix Zij can be related to the dressed charges by the algebraic relations: 

1 



^21 + Z21 



Zll = 62(fco)^12 

^22 ~ [2^22 (Ao)] ^ 

-ei2(fco)/62(Ao) 



(D.ll) 



while the remaining combinations are expressed in terms of ais and CT2s by 
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1 



■^12 — — :r 



21 



^21 ^ ^21 



2 

sin ko 
— sin fcn 



dt 



dXais{X) - 
1 



(iAcris(A) 



2n 



27r 



i^i(i-A)(T2.(A) 



(D.12) 



The two functions (Tis(A) and cr2s(A) satisfy the foUowing equations: 



ZTT 



if2(A-A')ai,(A') 

1 rl\' 
^2. (A) = - — i^2(A - Ao) - / ^if2(A - A')a2,(A') 
271- J-x„ 27r 



^° dA' 

^0 J\/ 



(D.13) 



By s ubstituting Eq. (D.7) into Eq. (D.6) and evaluating the quantities in brackets by means of Eq. (D.S) we get Eq. 
In conclusion, the relevant equations for the one hole problem in the Hubbard model are: 



— Xi(A-sinfc)cos2fc- / —K, 
-TT 27r J_^^ 27r 



(A-A')e,(A') 



— i4ri(sin fc 
^Ao 27r 



-A)e,(A) 



(D.14) 



which give the spin and charge excitation energies. Obviously, the holon energy is — —ec- The external fields 
and he are chosen in such a way that e(Ao) = and e(fco) = 0. The charge and spin distributions are given by: 



Pc{k) = 

The spinon and holon velocities are 



^-0Xi(A-sinfc)-y ^ ^K,{X-X')psiX') 
Ki{sink — X)ps{X) 



2tt 



cos k 



27rv, 



dX 



-A( 



2tt 



(D.15) 



1 dcsiX) 



2ttVc = - 



Ps(Ao) dX 
1 dec(A:) 



Pc{ko) dk 



(D.16) 



The cut-off Aq and the holon (dressed) momentum ko are related to the magnetization and to the position of the hole 
in the charge distribution by 



h 
L 

L 



dkpc{k) 



dXp,{X) 



(D.17) 



-An 



Finally, the elements of the dressed charg e mat rix (evaluated at the cutoff fcp and Aq) are given by Eqs. (D.IO) and 
the additional matrix is defined by Eqs. ( D.12 ) with the the help of Eqs. ( D.13 ). The finite size corrections to the 
energy and the correlation exponent Xg are then given by Eq. (|7.3|) in terms of the above defined quantities. 



2. Supersymmetric t-J Model 



The expressions for the t — J model are quite similar. In particular the formal size corrections (7.3) are the same 
while only the definition of the coefficients are different. Also the procedure closely follows that outlined for the 
Hubbard model. Therefore here we will just report the final expressions in Sutherland's representation |^,^ noting 
that in this formalism there is only one "charge" (i.e. the hole): Nc = 1- 
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The spin and charge excitation energies are given by the equations 



Uv) = hs- 29{2v) - I ^e{v - v')es{v') 



ec{w) = hc + 



dv . 
2^ 



27r 



(D.18) 



where hg and are chosen in such a way e(wo) — and e(wo) = 0. The kernel 9{x) = 2arctanx is defined according 
to Bares tt al |Q. Dot represents derivation with respect to the argument. The charge and spin distributions are 
given by: 



1 A, 



dv' ■ 



Psiv) ^ -9{2v) ~ —e{v-v')p4v') 

J-vo 271- 

r" dv ■ 

Pc{w)= / —2e{2w-2v)ps{v) 



and the spinon and holon velocities are: 



2'KV, 



2-K 



1 dts(v) 



(D.19) 



2'KVr = — 



PsM dv 
1 dtc{w 



Pciwo) dw 



(D.20) 



The cut-off vq and the holon rapidity wq are related to the magnetization and to the position of the hole in the charge 
distribution by 



h _ 
L " 

l + Ns 
L 

Other relevant quantities are the dressed charges: 

= 1- 
62 (w) = 



dvps{v)9{2wo — 2v) 



-VO 

■Vq 



dvps{v) 



"0 dv' ■ 

^e{v - V')^22{V') 



(D.21) 



'"° dv ■ 

— 29{2w - 2v)^22{v) 



(D.22) 



In the following, we will consider the dressed charges evaluated at the cut-off: ^22 = ^22{vq) and ^12 = ^12(^0). 
Finally, the elements of the Z-matrix are defined by: 



1 r'° dv 
Zii = -+ / —(7i{v)e{2wo ~ 2v) 

2 J-vo 

e{2vo - 2wo) - e{2vo + 2wo) H dv 



z 



1 



27r 



2ti 



(Ti{v)[e{vQ - v) - e{vfi + v)] 



^21 ^ ^21 — 

where the function ai {v) satisfies the following equation: 

1 • 



(Ti{v) ^ -e{2v - 2wo) - / —e{v -v')(Ti{v') 



dv' ■ 



27r 



(D.23) 



(D.24) 
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FIGURE CAPTIONS 



Fig. 1 Exact spectral function A{uj) (in units of the hopping t) for the single hole in the Ising model for two 
values of J. The heavy line represents the S contribution and its height is proportional to the quasiparticle 
weight. 

Fig. 2 Overlap square Z betwee n the Heisenberg ground state on a ^ = i — 1 site ring and the ground state 
of the effective spin Hamiltonian ( |3.7D corresponding to a single hole in the L site t — J model at J = At and 
p = n/2. Lanczos diagonalization has been performed on even chains with L < 26. The dashed line is a guide 
to the eye. 

Fig. 3 Overlap square Z b etw een the Heisenberg ground state on a i site ring and the ground state of the 



effective spin Hamiltonian (3/7) of a single hole in the L site t — J model a.t t — Q (static limit). Lanczos 
diagonalization has been performed on even chains with L < 26. The dashed line is a parabolic fit of Lanczos 
data. The full line shows the expected asymptotic slope of the curve on the basis of conformal field theory 
(2Xo = 3/8). 

Fig. 4 Leading exponent of the size dependence of the overlap C, in the t ~ Jxy model as a fun ction of the 
number N of up spins. Lines represent the num erical evaluation of the determinant of the matrix ( |6.11 ). Dots 



are the analytical prediction based on Eq. (6.7) 



Fig. 5 (a) Spinon term Z{Q) for the t— Jxy model computed as the determinant of the matrix (6.16) at r = 
for magnetization /i = (squares) and fi — 1/A (triangles). Two cases corresponding to about 400 and 800 
particles are shown by open and full symbols respectively, (b) Logarithmic plot of Z{Q) nea r th e singularity. 
The lines indicate the slope predicted on the basis of the finite size corrections to the energy (|6.7[). 



Fig. 6 Asymptotic behavior of the dynamic al q uasiparticle weight Z{R, r) in the t — Jxy model at /i = — 1/ 4 



calculated as the determinant of the matrix (|C/7|) in a system of 201 particles. Zo{R, t) is defined in Eq. (3.18) 
Data refer to more than a thousand different points in the (i?, r) plane. 

Fig. 7 Critical exponents Aq (full squares) and A_i (open squares) as a function of the momentum of the 
hole p for several magnetizations ji — —0.4, —0.3, —0.2, —0.1, 0.1, 0.2, 0.3, 0.4 from (a) to (h) respectively, in the 
Hubbard model aXU ^ At. 



Fig. 8 Loci of the singularities of the spectral function in the (p, w) plane (5.56) for the Hubbard model at 
U — At. Energies are measured from the bottom of the holon band. Two magnetizations are shown: /i = and 
jjL — 1/4. Full (dashed) lines correspond to critical exponents X smaller (larger) than 1/2. In the latter case, 
according to Eq. (|5.57|) , the divergence disappears. 



Fig. 9 Leading exponent A™ of the finite size scaling of the overlap C, defined in Eq. (5.36) as a function of 
the hole momentum p. Data refer to the t — J model at /i = ±1/4 and J — 2t in panel (a) and J — t in panel 
(b). Solid lines: analytical results obtained from the finite size corrections to the energy in the Bethe ansatz 
solution. Full (open) dots: Power law fit of Lanczos data for = 1/4 (^ = —1/4). Numerical data (full squares) 
obtained with the same fit for the non integrable point J = t, are shown for comparison only in the /i = 1/4 
case. The dashed line connecting these points is a guide to the eye. Diagonalizations have been performed on 
rings with L = 8, 16, 24, 32. 



Fig. 10 Breakdown of the completeness condition (8.3) in the single spinon approximation for the t — Jxy 



model at two magnetizations. L is the length of the chain. 

Fig. 11 Spinon function Z{Q) for the Hubbard model at /it = according to the single spinon approximation, 
(a) refers to U — At and p = pp', (b) to U — 8t and p = pp; (c) to U — At and hole momentum p — 0; (d) to 
U — 8t and p — 0. The hole Fermi momentum pp — tt/2 corresponds to the bottom of the (down) hole band. 
Z{Q) vanishes inside the spin up Fermi sphere. Symbols represent Lanczos data for different system sizes: L = 6 
(open circles), i = 10 (open squares), L = lA (full circles). 

Fig. 12 Same as Fig. 11 at magnetization fj, ~ 1/4 corresponding to pp — n/A. Symbols represent Lanczos 
data for different system sizes: L ~ 8 (open circles), i = 12 (open squares), i = 16 (full circles). 
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• Fig. 13 Hole spectral function in the Hubbard model at U = At as a function of frequency measured from the 
bottom of the holon band. Energy is in units of the hopping t and the calculation has been performed in single 
spinon approximation. Panel (a) corresponds to p = and ^ — 0, (b) to p = and — 1/4, (c) to p = pp and 
/X = 0, (d) to p = pf and /x = 1/4. Vertical lines show the location of singularities. Dashed lines identify the 
divergences induced by band structure effects; dotted lines show the frequencies of the non trivial singularities 
reported in Fig. 8. 



38 



T 1 r 



T 1 r 



T 1 r 



0.5 



0.4 



0.3 



0.2 



0.1 



— 
0.5 - 



0.4 



0.3 - 



0.2 



0.1 







I I I 



J=6t 



H — \ — h 



J=2t 



I I I 



I I I 



w/t 







Sorella and Parola Fig. 1 



"I 1 1 1 1 1 1 1 1 1 1 1 1 r 



0.95 



0.9 



0.85 - 



0.8 



J I I I I I I I I I I I I I 







0.05 0.1 
1/L 



0.15 



Sorella and Parola Fig. 




Sorella and Parola Fig. 3 




Sorella and Parola Fig. 4 




Sorella and Parola Fig. 



3 



1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 

(b) : 


■ ~~~~ - 

■ ^ 




a' ~ _ 

A 

1 1 1 1 1 


■ ■ — 

■ — 

* 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 



-6 -5 -4 -3 

In I Q/k, - 1 I 



Sorella and Parola Fig. 5b 



0.7 



"I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



0.65 



o 

isi 



0.6 



0.55 - 



A6=-l/4 



0.5 



I I I I I I I I I I I I I I I I I I I 
20 40 60 80 

R-VT 



1 



100 



Sorella and Parola Fig. 6 




Sorella and Parola Fig. 7 




Sorella and Parola Fig. 8 




Sorella and Parola Fig. 9 




Sorella and Parola Fig. 10 



isi 4 







1 1 1 1 1 1 1 

~ 1 □ 1 

~ 1 O 1 
• 


1 1 1 1 1 1 1 

(b) 

1* 1 ~ 

,□ 1 _ 
1 o 1 

1 1 1 T •^^ 1 1 


1 1 1 T ^ 1 1 1 

- ! • ! (c) 

1 ^ 1 
1 o 1 

1**1 
1* *i 

1 o, ,o 1 

1 1 1 1 1 1 1 


1 1 1 1 1 1 1 
(d) 1 1 

I* W 

In g □! 

1 O^n* *n^O 1 

1 1 1 1 1 1 1 



7t/2 7t 37t/2 7t/2 " 37T/2 2tt 
Q Q 



Sorella and Parola Fig. 







' 1 ' 1 ' 

(a) 

1 o 1 
1 * 1 

1 1 i r^^i 1 1 


' 1 ' 1 ' 

(b) 

!• 1 

1 o 1 - 
• 

1 1 1 1 


(0 

o o 

— 1 ° ° 1 

• • 




(d) 1 1 





7t/2 7t 3tt/2 7t/2 " 3tt/2 Stt 
Q Q 



Sorella and Parola Fig. 




Sorella and Parola Fig. 13 



